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This  report  presents  a  dynamic  gust  analysis  of  the  Boeing  Model  720B 
airplane  and  outlines  two  procedures  for  assessing  gust  design  strength 
for  future  civil  transports.  This  work  was  conducted  under  subcontract 
for  the  Lockheed-California  Company  in  support  of  their  study  contract 
with  the  Federal  Aviation  Agency  (FAA)  to  develop  a  power  spectral  gust 
design  procedure  for  civil  aircraft. 

The  procedures  outlined  are  based  on  two  approaches.  The  first  is  the 
design  envelope  approach,  the  second  is  the  flight  profile  approa  h.  In 
the  design  envelope  approach,  certain  flight  conditions  are  established 
by  successive  analyses  to  determine  the  most  critical  flight  condition 
for  a  given  portion  or  component  of  the  airplane  for  either  vertical  or 
lateral  gust  loads.  The  one-factor  level  flight  load  added  to  the  rms 
load  times  a  particular  constant  will  .lust  equal  the  limit  design  strength 
of  the  component.  This  constant,  <rv  i)<p  represents  an  effective  gust 
intensity  which  would  just  stress  the  structure  to  its  limit  design 
strength. 

The  object  of  the  flight  profile  approach  is  to  determine  the  expected 
number  of  flight  hours  that  the  airplane  could  be  operated  before  the 
limit  design  strength  of  any  of  its  major  components  would  be  exceeded. 

The  flight  profile  approach  requires  a  description  of  airplane  operation 
in  terms  of  flight  profiles  that  best  typify  the  airplane  usage.  A  sep¬ 
arate  power  spectral  analysis  is  conducted  for  each  of  the  profile  condi¬ 
tions.  In  addition,  a  description  of  the  atmosphere  applicable  to  the 
condition  altitude  is  determined.  From  this  information,  the  expected 
number  of  hours  required  to  exceed  the  limit  strength  is  computed. 

The  720B  airplane  was  studied  for  both  concepts,  first,  by  using  the 
bending  moment  on  the  wing,  fuselage,  and  vertical  tail  as  indicies  of 
their  strengths.  This  procedure  was  used  to  locate  critical  flight  con¬ 
ditions,  critical  portions  of  the  structure,  and  to  obtain  preliminary 
values  of  aw  by  the  design  envelope  approach,  and  expected  hours  to 
fly  to  exceed  limit  design  strength  by  the  flight  profile  approach. 

The  second  step  was  to  study  the  more  critical  structural  areas  using  a 
joint  probability  stress  analysis  approach  developed  by  Boeing.  The 
critical  values  of  0^  and  hours  to  fly  to  exceed  limit  design  strength 
of  each  structural  element  in  the  critical  area  were  determined.  The 
results  of  these  analyses  indicate  that  the  design  envelope  and  flight 
profile  approaches  give  results  that  are  in  general  agreemer*-. 
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Background 


Gust  load  formulas  embodying  a  number  of  simplifying  assumptions  have 
been  used  in  the  United  States  for  over  30  years  to  calculate  gust  loads 
for  airplanes.  The  underlying  concept  in  the  use  of  a  gust  formula  is 
that  measured  airplane  center  of  gravity  accelerations  due  to  continuous 
turbulence  can  be  used  to  derive  effective  or  derived  gust  velocities  for 
specified  gust  shapes.  These  derived  gust  velocities  for  a  specified 
gust  shape,  in  turn,  can  be  used  to  calculate  the  accelerations  on  other 
airplanes  by  reversing  the  process.  Therefore,  a  gust  load  formula  re¬ 
lates  the  peak  gust  induced  accelerations  on  a  given  airplane  to  the 
peak  accelerations  expected  on  other  more  or  less  similar  airplanes  for 
flight  through  the  same  continuous  rough  air.  It  is  clear  that  when  newer 
airplanes  differ  significantly  from  past  or  present  successful  airplanes 
on  which  acceleration  measurements  have  been  made,  the  simple  gust 
formula  approach  may  not  be  adequate. 

When  a  proposed  airplane  design  or  mode  of  operation  does  differ  signifi¬ 
cantly  from  past  practice,  the  aircraft  manufacturer  and  the  certifying 
agency  must  cope  with  the  problem  of  providing  adequate  strength  for  gust 
loading.  Ideally,  the  risk  of  exceeding  the  limit  design  gust  loads  per 
flight  hour  for  the  new  airplane  should  be  just  equal  to  the  risk  of 
exceeding  the  limit  design  gust  loads  per  flight  hour  for  the  older  proven 
airplanes. 

An  evaluation  to  establish  gust  loads  for  new  airplane  designs  should 
rely  on  rather  detailed  theoretical  dynamic  analyses  that  can  adequately 
describe  the  behavior  of  the  new  airplanes  relative  to  the  older  proven 
models.  A  generally  accepted  analysis  technique  that  is  used  to  account 
for  differences  in  airplane  dynamic  responses  and  in  modes  of  operation 
is  power  spectral  analysis.  The  theoretical  technique  for  using  gener¬ 
alized  harmonic  analysis  or  power  spectral  analysis  methods  in  the  air¬ 
plane  gust  loads  problem  have  been  developed  and  summarized  by  the 
National  Aeronautics  and  Space  Administration,  (NASA). 

The  contracting  agency,  Federal  Aviation  Agency,  (FAA),  has  assumed  that 
the  manufacturers  of  given  airplanes  can  best  analyze  their  own  designs, 
since  they  are  most  familiar  with  the  developmental  wind  tunnel  and 
structural  testing,  the  detailed  construction,  and  the  manner  in  which 
the  airplanes  have  been  used  in  day  by  day  service .  After  an  evaluation 


*Numbers  in  parentheses  refer  to  items  in  the  list  of  references. 
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of  detailed  proposals  submitted  for  this  study,  the  FAA  contracted  with 
the  Lockheed-California  Company,  Burbank,  California,  to  study  the 
Lockheed  Eiectra,  Model  l88,  a  medium  weight,  straight  wing,  4-engine, 
turboprop  airplane  and  the  Lockheed  Constellation,  Model  749,  a  medium 
weight,  straight-wing,  4-engine,  reciprocating  engine  airplane.  In 
addition,  the  FAA  arranged  to  have  The  Boeing  Company,  Airplane  Division, 
participate  in  the  study  under  a  subcontract  with  Lockheed  to  conduct  a 
study  of  the  Boeing  Model  720B,  a  medium  range,  swept  wing,  4-engine, 
jet  transport  airplane .  A  view  of  the  720B  airplane  in  flight  is  shown 
in  Fig.  1. 


Approach 

Even  though  the  basic  power  spectral  technique  is  well  defined,  the 
problem  of  establishing  specific  gust  design  levels  for  new  airplanes 
requires  that  a  concerted  effort  be  made  to  establish  the  relationship 
between  the  stresses  induced  in  critical  structural  elements  of  existing 
successful  airplanes  and  the  design  strength  of  these  elements.  In  this 
manner,  an  estimate  can  be  made  of  the  expected  number  of  flight  hours 
that  proven  airplanes  can  be  expected  to  operate  before  the  limit  design 
strengths  of  their  various  structural  elements  are  exceeded.  This  ap¬ 
proach  requires  an  intimate  knowledge  of  the  airplane  usage,  and  is 
commonly  referred  to  as  the  flight  profile  approach.  A  second  approach 
involves  computation  of  the  statistical  probability  of  exceeding  the 
limit  strengths  of  critical  structural  elements  when  the  airplane  _  being 
flown  in  a  manner  to  most  highly  stress  that  portion  of  the  airplane 
containing  the  critical  elements.  The  level  of  gust  intensity  at  a 
given  altitude  that  would  just  stress  the  structural  element  in  question 
to  its  limit  strength  provides  a  measure  of  gust  intensity  that  might 
logically  be  used  on  a  new  design  for  the  same  altitude.  This  second 
approach  is  referred  to  herein  as  the  design  envelope  approach.  Both  of 
these  approaches  are  exercised  in  the  analyses  that  follow. 

The  approach  of  establishing  gust  design  levels  for  new  airplanes  by 
strength  analyses  of  older  successful  airplanes  raises  the  question  of 
sensitivity  of  the  analysis  procedure  to  the  various  parameters  used  in 
the  analyses.  Therefore,  the  results  are  studied  to  determine  the  sen¬ 
sitivity  of  the  analyses  to  variations  in  aerodynamic  parameters,  the 
mathematical  model,  and  in  assumed  airplane  usage. 


NOMENCLATURE 


x ,  £  Shear  stress,  psi 

y  ,  f  Axial  stress,  psi 

xo,  £o>  Vo>  *o  Steady  stress  values,  shear  and  axial  stress, 

respectively,  psi 

Allowable  shear  stress,  psi 

Ft  Allowable  axial  stress,  psi 

Z  Stress  vector, 

Z  Stress  velocity  vector, 

oX)  Oy>  ° a  )  °ft  Root  mean  square  values  (standard  deviations) 

for  x,  y,  a  s  and  >  respectively 

$x(w)>  $y(  0) )  Power  spectral  density  functions  for  random 

process  x(t)  and  y(t),  respectively 

0Xv(  w )  Cross  power  spectral  density  function  for 

x(t)  and  y(t) 


0) 

P 

a  ,/3 

p  (x),  p(y) 

P  (x,y) 

f  (x,  a  ,y>  j8) 
MS 

p  (ms< o) 


Circular  frequency,  rad.  per  sec. 

Correlation  coefficient  for  x(t)  and  y(t) 
Time  rate  of  change  of  x  and  y,  respectively 
Probability  densities  of  x  and  y 
Joint  probability  density  of  x  and  y 
Probability  density  of  x,  a  ,y,  ft 
Margin  of  safety 

Probability  that  MS  is  less  than  zero  or 
percent  time  that  MS  is  less  thsn  zero 
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P(MS>0) 


C 

N 

[] 

LJ 

Nc 


av>  °v 

$(*v> 

V 

G 


No 


Probability  that  MS  is  greater  than  or 
equal  to  zero 

A  curve  on  xy-plane 

Unit  vector  normal  on  C 

Column  matrix 

Row  matrix 

Total  number  of  passages  across  an  arbitrary 
curve  per  second  or  per  foot  traveled 

Root  mean  square  gust  velocity,  ft.  per  sec. 

Probability  density  distribution  of  av 

Velocity,  ft.  per  sec. 

Expected  exceedances  of  limit  design  strength 
per  hour 

Root  mean  square  stress  response  for  ctw  of 
unity,  psi  per  ft.  per  sec. 

Humber  of  times  per  unit  time  or  distance 
that  a  time-history  crosses  its  mean  value 
with  positive  slope;  or  negative  slope 


DESCRIPTION  OF  THE  ATMOSPHERE 


Gust  Power  Spectrum 

Certain  statistical  estimates  for  Loads,  accelerations  or  stresses  for 
flight  vehicles  can  be  obtained  by  generalized  harmonic  analysis  tech¬ 
niques  if  sufficient  statistical  information  is  available  for  the  random 
forcing  environment.  If  it  can  be  assumed,  as  has  been  done  fir  the 
current  study,  that  (l)  the  turbulence  is  essentially  frozen  and  the  air¬ 
plane  penetrates  or  passes  over  the  turbulence  much  as  an  automobile 
would  travel  a  rough  road,  and  (2)  if  it  is  assumed  that  the  spanwise 
effects  of  the  turbulence  are  not  of  major  importance,  then  a  one¬ 
dimensional  power  spectral  approach  is  adequate. 

If  it  is  further  assumed  that  the  turbulence  is  isotropic;  that  is, 
statistically  invariant  for  any  rotation  of  translation  of  the  coordinate 
axis  system;  then,  a  single  power  spectral  density  function  can  be  used 
for  separate  vertical  and  lateral  analyses  for  a  given  flight  condition. 
The  assumption  that  there  is  no  significant  aerodynamic  or  inertia 
coupling  between  vertical  and  lateral  airplane  degrees  of  freedom  is 
made  when  it  is  assumed  that  separate  or  uncoupled  vertical  and  lateral 
dynamic  analyses  will  suffice.  All  of  these  rather  restrictive  assump¬ 
tions  have  been  made  for  the  analyses  discussed  herein;  and  a  constant 
parameter  normalized  power  spectral  density  function  was  chosen  to  rep¬ 
resent  atmospheric  turbulence  for  all  altitudes. 

A  meeting  was  held  with  NASA  and  FAA  personnel  at  Langley  Field  on 
March  10,  196U,  to  discuss  the  atmospheric  turbulence  description  to  be 
used  in  this  program. (3)  These  discussions  resulted  in  the  choice  of 
an  isotropic  turbulence  power  spectral  density  function  and  NASA's  con¬ 
sent  to* provide  operational  VGH  data  for  the  three  airplane  types  being 
studied.  The  VGH  data  was  provided  to  aid  Lockheed  in  determining  a 
probability  distribution  for  root-mean-square  (rms)  gust  velocity  to  be 
used  with  the  selected  spectrum  shape.  These  two  functions  provided  a 
sufficient  atmospheric  turbulence  model  for  the  overall  study. 
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The  atmospheric  turbulence  of  spectrum  chosen  was  as  follows: 
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a  2L 
w 


n 


1  +  |  (1.339  LI))2 

i  +  (1.339  Lfi)2* 


(1) 


A  scale  of  turbulence,  L,  of  2,500  feet  was  selected.  This  spectrum 
shape,  designated  as  the  "isotropic  turbulence,"  Case  I  spectrum,  has 
been  extensively  evaluated  by  NASA  and  appears  to  provide  a  good  fit  to 
their  experimentally  determined  gust  spectra. (2)  The  2,500  ft.  scale  of 
turbulence  is  most  probably  low,  and  would  be  expected  to  give  a  somewhat 
conservative  or  high  number  of  load  or  stress  exceedance  counts  per  unit 
time.  NASA's  experience  indicates  a  scale  of  turbulence  for  the  above 
power  spectrum  in  the  order  of  3>000  to  6,000  ft.  A  comparison  of  the 
Case  I  spectrum  shape  chosen  for  this  study,  and  the  commonly  used  Case 
II  isotropic  turbulence  spectrum  with  a  scale  of  turbulence  of  1,000  ft. 
is  shown  in  Fig.  2. 


Distribution  of  RMS  Gust  Velocities 


The  power  spectral  analysis  approach  used  herein  is  applicable  only  to 
stationary  Gaussian  continuous  turbulence,  but  atmospheric  turbulence  is 
not  statistically  stationary  or  Gaussian  over  long  distances.  The 
statistical  quantities  used  to  describe  turbulence  vary  with  altitude, 
wind  condition,  terrain  roughness,  temperature  gradient,  season  of  the 
year,  and  a  host  of  other  variables.  However,  it  has  been  observed  that 
the  power  spectrum  shape  from  1,000  to  1*0,000  ft.  above  the  terrain  is 
reasonably  invariant.  As  a  result,  NASA  has  proposed  that  atmospheric 
turbulence  be  considered  locally  Gaussian  and  stationary  and  that  the 
total  experience  of  flying  through  rough  air  be  considered  to  be  made  up 
of  a  range  of  exposures  to  turbulence  of  various  intensities  all  using 
the  same  shape  power  spectrum.  Thus,  they  have  proposed  that  a  statis¬ 
tical  distribution  of  rms  gust  intensities  be  considered. (^)  The  pro¬ 
posed  probability  density  function  for  rms  gust  velocity  is, 


A,  s 


%  4b, 


(2) 


Eq.  (2)  leads  to  an  expression  for  the  number  of  times  per  unit  time  that 
a  given  load  level,  y,  is  exceeded, 


G(y)  =  PJe'y/¥  +PN  e~y/b2^ 
-i*  o  2  o 


(3) 


where, 

N0  number  of  times  per  unit  time  that  the 

y-time  history  crosses  its  mean  value 
with  positive  (or  negative)  slope 

P  ,  P  fraction  of  time  or  distance  flown  in 

"  '  nonstorm  (light)  turbulence,  and  fraction 

of  time  or  distance  flown  in  storm  (heavy) 
turbulence,  respectively 

b^  b2  turbulence  scale  parameters  for  nonstorm 
and  storm  turbulence,  respectively 

A  rms  value  for  the  load  quantity,  y,  per 

unit  rms  gust  velocity 

The  zero  crossings  per  unit  time,  N0,  is  a  function  of  the  shape  of  the 
gust  power  spectral  density  function.  As  a  result,  when  a  new  spectrum 
shape  is  chosen,  the  values  for  P-j_,  P2,  b^,  and  b2  must  be  re-evaluated 
to  provide  load  results  which  agree  with  flight  data 

Lockheed  investigated  the  turbulence  parameters  Pi,  P2,  bi,  b2  and  ar¬ 
rived  at  the  relationships  shown  in  Figs.  3  and.  4.  These  turbulence 
parameters  were  determined  not  only  from  a  substantial  sample  of  airline 
VGH  data,  but  also  from  data  obtained  from  military  operations  at  higher 
altitudes.  Most  of  the  basic  data  was  obtained  from  a  summary  report 
prepared  by  the  Douglas  Aircraft  Company,  Inc.,  Long  Beach,  California, 
under  contract  with  the  Flight  Dynamics  Laboratory,  Wright-Patterson  Air 
Force  Base,  Ohio.  W 
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PROPORTION  OF  FLIGHT  TIME  OR  DISTANCE 

FIGURE  3.  PROPORTION  OF  FLIGHT  TIME  OR 
DISTANCE  IN  STORM  AND  NONSTORM  TURBULENCE 
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FIGURE  4.  TURBULENCE  SCALE  PARAMETERS 
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ANALYSIS  CONDITIONS 


Estimated  Airline  Usage  of  72QB  Airplane 

The  airplane  configurations  analyzed  during  this  study  were  selected  to 
provide  statistical  loads  and  stress  information  for  use  with  two  dif¬ 
ferent  design  approaches.  One  approach,  flight  profile  analysis,  re¬ 
quires  information  based  on  actual  airplane  usage  and  reflects  the  speeds, 
altitudes,  and  gross  weights  experienced  during  routine  airline  operation. 
The  second  approach,  design  envelope  analysis,  requires  the  generation 
of  statistical  loads  estimates  for  critical  design  points  on  or  within 
the  speed-altitude  design  envelope . 

A  survey  of  the  airline  usage  of  720  and  720B  airplanes  was  made  to 
determine  the  airplane  usage  most  appropriate  for  flight  profile  analy¬ 
sis.  Three  airlines,  whose  route  structures  were  representative  of  the 
fleet  usage,  were  chosen  for  the  survey.  These  airlines  operate  a  total 
of  thirty-seven  720/720B  airplanes  or  30  percent  of  the  entire  720/720B 
fleet.  The  survey  consisted  of  compiling  the  usage  of  these  airplanes 
in  terms  of  the  total  profiles  flown  for  particular  profile  lengths. 

The  profile  intervals  chosen  varied  from  100  to  2,600  n.  miles  in  100  n. 
mile  intervals.  One  flight  profile  or  segment  is  comprised  of  a  single 
takeoff,  climb-out,  cruise,  and  descent.  The  results  of  this  survey  are 
presented  in  Fig.  5(a)  in  terms  of  profile  length  as  percent  of  the  total 
profiles  flown.  A  comparison  of  these  results  with  the  results  of  the 
NASA  VGH  data  for  707-300  airplanes  (3)  is  shown  in  Fig.  5(b) .  The 
comparison  is  made  in  terms  of  the  percent  of  total  profiles  flown  versus 
profile  duration  in  minutes.  It  can  be  seen  that  there  is  a  definite 
similarity  between  the  operation  of  707-300  airplanes  and  720/720B.  air¬ 
planes  for  the  shorter  flights.  As  would  be  expected,  the  longer  range 
707-300  shows  a  greater  usage  over  longer  segments  than  720/720B 
airplanes. 

The  selection  of  representative  analysis  profiles  was  based  on  a  distri¬ 
bution  of  time  spent  on  the  various  flights.  .  This  distribution,  shown 
in  Fig.  5(c),  differs  from  the  percent  of  segments  flown,  Fig.  5(a), 
because  each  segment  was  weighted  by  the  segment  flight  time.  The  flight 
segment  duration  information  in  Fig.  5(c)  was  used  to  establish  the  re¬ 
presentative  flight  profiles  for  the  analysis  of  the  720B  airplane. 

Values  of  indicated  airspeed  versus  altitude  based  on  the  NASA  VGH 
data  (3)  were  used  for  the  flight  profile  studies.  A  comparison  of  the 
recommended  airspeeds  with  the  average  operational  speeds  is  presented  in 
Fig.  6.  This  comparison  shows  that  the  actual  operating  airspeeds  are 
less  than  the  recommended  speeds  at  the  lower  altitudes,  while  at  the 
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higher  altitudes  there  is  better  agreement  between  the  recommended  and 
operational  speeds.  The  statistical  values  for  indicated  airspeed  were 
used  in  the  flight  profile  analyses  in  order  to  be  as  consistent  as 
possible  with  the  operational  experience. 


Flight  Profiles  for  720B  Analyses 

Five  profiles  were  chosen  from  Fig.  5(c)  for  the  dynamic  analyses.  The 
five  basic  and  alternate  profiles  are  shown  in  Figs.  7  through  11.  Pro¬ 
files  I  through  V  have  ranges  of  350.,  600,  1000,  1600,  and  2500  n.  miles 
for  flight  durations  of  55,  85,  135,  205  and  325  minutes,  respectively. 

The  variations  of  gross  weight  and  speed  are  also  given  in  Figs.  7  through 
11.  The  gross  weight  at  the  beginning  of  each  flight  profile  is  based 
on  the  prescribed  fuel  loading  necessary  to  accomplish  the  flight  plus 
10,000  lb.  of  reserve  fuel  for  holding  or  alternate  field  requirements. 

The  cargo  compartments  were  considered  to  be  fully  loaded  on  a  space 
limited  basis,  and  the  passenger  load  factor  was  assumed  to  be  55  percent. 
This  load  factor  was  the  1963  average  for  720  and  720B  operations  by 
United  States  airlines. 


Flight  Profile  Analysis  Conditions 


A  summary  of  the  flight  profile  analysis  conditions  is  given  in  Tables  1 
and  2.  Conditions  1,  2,  and  3  were  selected  to  represent  the  climb 
portion  of  Profile  III,  Fig.  9.  Results  from  analysis  for  these  condi¬ 
tions  allowed  a  study  to  be  made  of  the  variations  of  load  and  stress 
exceedance  data  for  climb  at  an  average  climb  gross  weight.  Parametric 
variations  in  gross  weight  at  the  15,000  ft.  altitude,  Conditions  4,  5, 
and  6,  provided  information  to  account  for  the  climb  portions  of  other 
profiles.  Three  cruise  conditions,  Conditions  7,  8,  and  9,  were  selected 
to  examine  the  variations  in  loads  with  cruise  altitude  at  a  constant 
gross  weight.  These  points  are  indicated  on  Profile  II,  Fig.  8,  as 
Options  A  and  B,  and  on  Profile  V,  Fig.  11.  Two  additional  conditions, 
Conditions  10  and  11,  were  studied  for  the  same  35,000  ft.  altitude  at 
different  gross  weights  to  account  for  the  fuel  burn-off  during  cruise . 
Conditions  12  and  13,  for  altitudes  of  10,000  and  20,000  ft.,  respec¬ 
tively,  were  chosen  to  evaluate  the  descent.  The  descent  portions  of 
all  the  profiles  are  the  same;  therefore,  the  results  obtained  for  one 
descent  were  applied  to  all 'profiles .  , 
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FIGURE  8b.  FLIGHT  PROFILE  IIB 
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FLIGHT  PROFILE  ANALYSIS  CONDITIONS 
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Design  Envelope  Analysis  Conditions 


The  design  envelope  analyses  conditions  selected  initially  were  based  on 
conditions  used  in  the  airplane  design  loads  analysis  (5),  a  summary 
of  all  design  envelope  conditions  analyzed  for  the  vertical  and  lateral 
analyses  are  given  in  Tablet)  3  anti  4,  respectively.  The  first  design 
envelope  vertical  analyses  were  conducted  for  the  airplane  cruise  speed, 
Vq,  of  375  kt.  equivalent,  Mach  .873;  at  an  altitude  of  22,000  ft.  Four 
gross  weight  configurations  were  analyzed  for  this  speed  and  altitude  to 
determine  the  most  critical  condition.  The  conditions  were  operating 
weight  empty,  119,700  lb.:  maximum  zero  fuel  weight,  156,000  lb.;  fuel 
transfer  weight,  186,000  lb.;  and  maximum  flight  weight,  221,600  lb.  The 
conditions  that  were  most  critical  were  then  examined  at  Vq  speed,  375 
kt.  at  5,000  and  15,000  ft.  and  at  Mq  speed,  Mach  No.  .90  at  30,000  and 
40,000  ft.  to  obtain  the  load  and  stress  variations  with  changes  in 
altitude .  At  the  critical  altitudes  determined  from  the  foregoing  result 
three  speed  variations  were  conducted.  The  speeds  chosen  for  the  varia¬ 
tions  were  VB,  253  kt.  EAS;  300  kt.  EAS  and  VD,  445  kt.  EAS  to  Mach  .95. 
Combining  the  results  from  the  speed,  altitude  and  weight  configuration 
variations,  the  boundary  of  minimum  margins  was  determined.  A  tabulation 
of  the  design  envelope  analyses  conditions  for  the  vertical  analysis  is 
given  in  Table  3.  These  conditions  are  also  indicated  on  the  speed 
altitude  diagram  in  Fig.  12. 

The  design  envelope  conditions  selected  to  isolate  the  most  critical 
point  for  the  lateral  analyses  were  slightly  different  from  those  for  the 
vertical  analyses  because  the  lateral  response  is  not  sensitive  to  the 
specific  details  of  fuel  loading.  A  high  gross  weight,  which  tends  to 
give  lower  Dutch  roll  stability,  was  selected  and  six  solutions  were 
obtained  at  Ve  for  altitudes  ranging  from  15,000  to  43,000  ft. 

Additional  conditions  were  analyzed  to  evaluate  the  effects  of  other 
parameter  variations  such  as  wing  bending  and  torsional  stiffness,  fuel 
loading,  number  and  limits  of  discrete  frequencies  used  in  the  analysis, 
yaw  damper  operation  and  lateral  aerodynamic  stability  parameters.  A 
summary  to  the  parametric  variation  conditions  is  given  in  Table  5» 


DESIGN  ENVELOPE  CONDITIONS  -  VERTICAL  ANALYSES 
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ANALYSIS 


General  Remarks 


The  Boeing  Mode]  720B  airplane  was  used  to  idealize  the  basic  mathemat¬ 
ical  models  for  all  analyses.  Fig.  1  shows  the  airplane  in  flight. 

This  airplane  requires  a  crew  of  three,  pilot,  copilot  and  flight  engi¬ 
neer,  and  has  a  maximum  range  of  3>300  miles.  The  airplare  normally 
cruises  from  15,000  to  ^0,000  ft.  at  gross  weights  in  excess  of'  200, 
000  lb.  and  at  speeds  up  to  615  miles  per  hour.  The  airplane  has  a 
wing  span  of  approximately  131  ft.,  and  an  overall  length  of  approxi¬ 
mately  137  ft.  A  more  detailed  3-view  diagram  of  the  airplane  is 
shown  in  Fig.  13. 

All  major  parts. of  the  airplane  except  the  horizontal  stabilizer  were 
considered  to  be  elastic  in  the  analyses.  Therefore,  a  rather  compre¬ 
hensive  mass  and  stiffness  description  of  the  airplane  was  required. 
Simple  beam  bending  theory  was  used  to  represent  the  stiffness  char¬ 
acteristics  of  the  major  components  such  as  the  wing,  forward  and  aft 
fuselage,  and  the  vertical  tail.  The  elastic  axes  were  generally  lo¬ 
cated  on  the  locus  of  shear  centers  along  each  component,  but  were  ad¬ 
justed  to  agree  with  information  obtained  in  previous  airplane  static 
testing. 

Structural  damping  for  the  various  airplane  flexible  modes  of  vibration 
was  taken  directly  from  whole  airplane  shake  test  results^1?'. 

Considerable  airplane  test  data  was  available  to  formulate  the  mathe¬ 
matical  models  that  were  used  in  the  analyses  of  the  720B.  The  method 
used  to  represent  the  detailed  aerodynamics,  mass,  and  stiffness 
characteristics  are  not  complex,  but  they  do  provide  a  solid  basis  for 
understanding  the  airplane's  dynamic  response  characteristics.  The 
airplane  responses  obtained  herein  by  analysis  show  good  agreement  with 
static  aeroelastic  results  for  slow  maneuvers  and  with  flutter  analyses 
for  the  higher  frequency  elastic  responses. 

Airplane  Mass  Data 


The  weight  and  inertia  data  for  the  analyses  were  originally  calculated 
from  released  detailed  drawings  and  have  been  updated  as  drawing  changes 
or  actual  weights  became  available.  The  airplane  fuselage  was  divided 
into  eighteen  panels  and  the  mass  properties  in  the  form  of  weight, 
center  of  gravity,  and  mass  moments  of  inertia  for  pitch,  roll  and  yaw 
were  determined  for  each  section.  Mass  data  for  the  operating  weight 
empty  (OWE)  configuration  and  for  full  and  partial  payloads  were  cal¬ 
culated  to  facilitate  selection  of  partial  payload  conditions. 


Um'iimJ  ►  DIMENSIONS  SHOWN  ARE  BASED 

UPON  GROSS  WEIGHT  OF  110,000 
POUNDS  AND  CENTER  OF  GRAVITY 
AT  15%  OF  MAC 


FIGURE  13.  3  -  VIEW  DIAGRAM  OF  720B 
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The  wing  structure  and  contents  with  the  exception  of  fuel  were  divided 
into  ten  spanwise  panels.  Each  of  these  panels  was  further  divided  into 
five  zones:  the  leading  edge,  front  spar,  inner  spar,  rear  spar  and 
trailing  edge.  The  weight  and  center  of  gravity  were  calculated  for  each 
zone  and  summed  to  give  the  total  panel  weight  and  center  of  gravity. 

The  total  panel  mass  moments  of  inertia  were  computed  by  rotation  and 
transfer  of  zone  results  into  axes  parallel  and  perpendicular  to  the 
wing  elastic  axis. 

The  vertical  tail  structure  and  contents  were  divided  into  seven  span- 
wise  panels  and  the  section  mass  properties  were  calculated  in  the  same 
manner  as  the  wing.  The  horizontal  stabilizer  was  considered  as  a 
single  panel  for  these  analyses.  The  lumped  mass  properties  for  the 
horizontal  stabilizer  were  obtained  by  combining  previously  panelled 
data  which  were  obtained  in  the  same  manner  as  the  wing  and  vertical 
tail  data. 

The  mass  properties  for  each  JT3D-1  engine,  nacelle,  and  nacelle  strut 
were  combined,  and  a  lumped  center  of  gravity  location  determined.  Then, 
the  nacelle  mass  moments  of  inertia  were  determined  for  axes  perpen¬ 
dicular  and  parallel  to  the  airplane  reference  axes. 

The  fuel  mass  properties  were  obtained  for  conditions  representing  fuel 
loadings  of  25,  50,  75  and  100  percent  of  full  fuel.  The  particular 
percentage  used  was  determined  by  the  flight  condition  being  analyzed. 
Fuel  weights  and  eg  positions  were  calculated  by  a  computer  program 
which  accounted  for  the  geometry  of  each  tank. 

The  accuracy  of  the  mass  properties  used  in  these  analyses  is  estimated 
to  be  within  1.0  percent  for  the  total  airplane  weight,  2.0  in.  for 
total  airplane  center  of  gravity  location,  and  3-0  percent  for  the  total 
airplane  mass  moments  of  inertia.  Also,  it  is  estimated  that  the 
panelled  mass  data  and  the  major  airplane  component  mass  data  are  ac¬ 
curate  to  similar  limit. 

Component  Stiffnesses 


The  stiffness  characteristics  of  each  major  component  of  the  airplane 
was  described  by  a  distribution  of  bending  stiffness,  El,  and  torsional 
stiffness,  GJ,  along  the  elastic  axis.  The  wing  section  properties 
were  computed  using  front  and  rear  spar  areas  and  all  inspar  skin  for 
both  the  upper  and  lower  surfaces.  Values  for  the  modulus. and  shear 
modulus  of  elasticity,  E  and  G,  were  10. 3(l0)°  and  3.8(lO)°  psi,  re¬ 
spectively.  The  body  section  properties  were  computed  using  stiffeners 
with  full  skin  effective  in  tension  and  2W  widths  of  skin  effective  in 
compression.  The  body  cutout  sections  were  handled  on  an  individual 
basis  by  special  analysis. 
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The  section  properties  for  determining  the  bending  stiffness  of  the  ver¬ 
tical  tail  were  computed  using  front  and  rear  spar  chord  area  with  2W 
widths  of  skin  assumed  effective  on  the  compression  side  of  the  rear 
spar  chord  in  trailing  edge  beam  areas.  Fully  effective  tension  skin 
was  assumed  between  the  front  and  rear  spars  and  for  balance  panel 
covers.  Near  the  body  intersection,  the  inspar  skin  and  front  spar 
areas  are  assumed  to  be  partially  effective  because  of  shear  lag.  The 
section  properties  for  the  determination  of  the  torsional  stiffness  of 
the  vertical  tail  were  computed  using  front  and  rear  spar  areas  and  all 
inspar  skin  for  both  surfaces. 


Aerodynamics  -  Vertical  Analyses 

The  stresses  induced  into  the  airplane  structure  for  unaccelerated  flight 
in  turbulent  air  result  from  the  sum  of  the  one  factor  aerodynamic  and 
inertia  forces  on  the  airplane,  the  gust  forces,  and  the  dynamic  re¬ 
sponse  forces.  The  one  factor  flight  loads  and  stresses  were  added  to 
the  respective  incremental  dynamic  loads  and  stresses  to  assess  the 
strength  capabilities  of  the  720B  airplane  for  flight  in  turbulent  air. 

The  one  factor  level  flight  loads,  that  is  bending, moments,  shears,  and 
torsions,  were  determined  by  aeroelastic  analyses.'1  All  of  the 

basic  aerodynamic  data  required  for  these  analyses  were  obtained  from  a 
series  of  wind  tunnel  tests  performed  for  the  entire  Mach  number  range. 
Wind  tunnel  pressure  model  test  results  were  used  to  establish  wing, 
fuselage,  an!  vertical  tail  airload  distributions.  The  section  lift  and 
moment  coefficients  were  corrected  for  model  flexibility  before  they 
were  used  for  full  scale  airplane  analysis,  and  were  later  refined  to 
obtain  final  agreement  between  the  aeroelastic  analyses  and  actual  air¬ 
plane  flight  load  survey  measurements. 

The  rigid  airplane  aerodynamic  section  coefficients,  corrected  to  pro¬ 
vide  agreement  with  flight  measurements,  can  be  used  to  obtain  an 
aeroelastic  solution  for  any  longitudinal  maneuvered  flight  condition. 
Simultaneous  equations  are  solved  to  obtain  an  aeroelastic  solution. 

These  equations  define  the  following:  (l)  wing  lift  distribution  on  10 
aerodynamic  sections  on  the  wing  semi-span,  (2)  airplane  lift  balance, 
and  (3)  the  airplane  pitching  moment  balance.  Specifically,  the  solution 
gives  the  elastic  wing  lift  distribution,  the  airplane  wing  root  angle 
of  attack,  and  the  balancing  tail  load. 

The  unsteady  aerodynamics  used  in  the  vertical  dynamic  analyses  are 
based  on  modified  two-dimensional  strip  theory. (”)  The  strip  theory, 
originally  developed  by  Theodorsen(9)  and  Kussner, (10/  has  been  modified 
to  include  aerodynamic  induction  effects(^).  These  induction  effects 
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account  for  the  aerodynamic  pressure  carry-over  between  wing  panels  and 
between  the  wing  and  horizontal  tail.  This  is  accomplished  by  using  a 
downwash  matrix  based  on  lifting  line  theory.  The  dynamic  downwash 
matrix  includes  pressure  carry-over  and  pressure  transmittal  functions 
to  provide  the  proper  magnitude  and  phasing  of  tne  carry-over  pressures. 
The  section  (strip)  aerodynamics  for  zero  frequency  are  made  to  agree 
with  the  comparable  aeroelastic  solution.  The  results  of  this  modifi¬ 
cation  have  been  compared  to  experimental  results  and  to  theoretical 
results  based  on  lifting  surface  theories;  satisfactory  correlation  was 
obtained. 

The  unsteady  aerodynamic  expressions  for  the  vertical  analyses  were  cal¬ 
culated  for  a  range  of  reduced  frequencies  to  adequately  define  the 
airplane  dynamic  responses.  The  gust  or  excitation  airforces  for  a 
sinusoidal  gust  field  were  also  expressed  in  a  modified  two-dimensional 
form  for  the  same  reduced  frequency  values.  The  aerodynamic  forces 
acting  on  the  fuselage  were  obtained  using  simple  quasi-steady  aerody¬ 
namic  expressions.  The  fuselage  was  first  divided  into  thirteen  stream- 
wise  panels,  seven  on  the  forebody  and  six  on  the  aftbody.  A  static 
pressure  distribution  was  then  used  to  determine  the  individual  panel 
lift  and  moment  coefficients.  The  panel  angles  of  attack  were  assumed 
to  arise  from  the  deflection  slopes  and  displacement  velocities  of  the 
fuselage. 

Aerodynamics  -  Lateral  Analyses 


The  lateral  dynamic  response  of  large' swept  wing  subsonic  jet  transports 
is  generally  dominated  by  the  low  damped  Dutch  roll  mode.  The  damping 
of  this  mode,  as  determined  by  dynamic  analysis,  is  very  sensitive  to 
the  estimated  ro . .  /y  stability  derivatives  used  in  the  equations  of 
motion.  Large  changes  in  the  lateral  gust  loads  result. from  small 
changes  in  particular  rotary  derivatives  for  flight  conditions  in  which 
the  Dutch  roll  damping  is  low. 

All  of  the  airforces  and  inertia  forces  were  referenced  to  the  airplane 
stability  axes.  The  stability  axes  were  determined  as  the  orthogonal 
set  attached  to  the  airplane  with  the  initial  orientation  of  the  x-axis 
in  the  wind  direction  for  undisturbed  and  unaccelerated  flight.  The 
airplane  angle  of  attack  for  one  factor  level  flight  was  determined  by 
a  separate  aeroelastic  analysis  for  each  condition. 

The  lateral  stability  derivatives  were  decomposed  into  wing,  fuselage, 
and  vertical  tail  component's.  The  pressure  distributions  for  the  fuse¬ 
lage  and  the  vertical  tail  were  determined  from  wind  tunnel  pressure 
model  test  data.  The  aerodynamic  side  force  and  moment  from  the  wing 
were  accounted  for  by  a  concentrated  side  force  and  a  moment  applied  to 


the  airplane  center  of  gravity.  The  integrated  force  and  moment  resulting 
from  the  aerodynamic  pressures  on  the  fuselage  and  vertical  tail  together 
with  the  concentrated  force  and  moment  from  the  wing  were  made  to  agree 
with  the  lateral  whole  airplane  derivatives. 

Generalized  Coordinates 


Each  dynamic  analysis  was  conducted  for  a  mathematical  model  which  was  in 
equilibrium  at  one  factor  level  flight.  When  such  a  system  is  subjected 
to  an  external  sinusoidal  gust  field,  it  responds  about  its  static  equi¬ 
librium  configuration.  This  gives  rise  to  lisplacements,  velocities  and 
accelerations  necessary  to  preserve  the  energy  balance  with  the  external 
gust  loading.  The  structural  deformation  of  the  airplane  can  be  de¬ 
scribed  as  a  sum  of  incremental  displacements  in  the  system's  natural 
modes  of  vibration.  These  deflected  shapes  and  whole  airplane  rigid 
body  displacements  are  used  as  the  so-called  generalized  coordinates  to 
describe  the  airplane  motions.  The  development  of  the  equations  of 
motion  by  the  Lagrangian  approach  using  these  generalized  coordinates  is 
discussed  in  Appendix  A. 

The  natural  modes  of  vibration  of  a  flexible  system  can  take  various 
forms  depending  upon  the  manner  in  which  the  system  is  restrained.  When 
the  airplane  is  in  flight  there  are  no  external  restraints,  and  the 
airplane  can  displace  with  complete  freedom.  There  are  several  methods 
by  which  this  so-called  free-free  state  can  be  described.  One  method  is 
to  define  the  cantilevered  vibration  modes  of  each  major  component,  such 
as  the  wing,  body,  and  nacelles  by  separate  orthogonal  mode  sets.  A 
limited  number  of  these  cantilevered  modes  along  with  the  rigid  airplane 
freedoms  are  used  as  the  generalized  coordinates  to  describe  the  overall 
airplane  response. 

A  better  set  of  generalized  coordinates  is  obtained  when  a  large  number 
of  cantilevered  modes  are  coupled  with  the  rigid  airplane  freedoms  to 
produce  a  set  of  unrestrained  free-free  airplane  normal  modes.  The  ad¬ 
vantages  of  this  approach  is  that  the  overall  airplane  response  can  be 
more  accurately  described  by  fewer  degrees  of  freedom  or  generalized 
coordinates  than  in  the  cantilevered  mode  approach.  The  free-free  normal 
mode  approach  was  used  for  the  720B  vertical  analyses,  and  the  canti¬ 
levered  mode  approach  was  used  for  the  lateral  analyses.  A  representative 
summary  of  cantilevered  modes  which  were  cou.  i.ed  with  airplane  vertical 
translation  and  pitch  freedoms  for  the  vertical  analyses  is  given  in 
Table  6.  Corresponding  representative  airplane*  free-free  modal  data  is 
given  in  Table  7*  A  summary  of _ cantilevered  modes  used  in  the  lateral 
analyses  is  given  in  Table  8. 


TABLE  6 


REPRESENTATIVE  CANTILEVERED  MODES 
VERTICAL  ANALYSES 


No.  of 

Cantilevered 

Modal  Frequency, 

cps 

Airplane 

Discrete 

Mode 

OWE 

Half 

Full 

Component 

Masses** 

Description 

Fuel 

Fuel 

Wing 

10 

1st  Bending 

1.31 

1.27 

.93 

2nd  Bending 

3.77 

3.51 

2.81 

3rd  Bending 

7.75 

6.59 

5.83 

4th  Bending 

10.67 

10.62 

9.28 

1st  Torsion 

2.52 

2.52 

2. 52 

2nd  Torsion 

5.19 

5.18 

5.19 

3rd  Torsion 

19.79 

19.28 

18.42 

Forward 

Fuselage 

7 

1st  Bending 

3-99 

3.71 

3.71 

Aft 

13* 

1st  Bending 

2.57 

2. 46 

2.46 

Fuselage 

2nd  Bending 

9.87 

8.62 

8.62 

Inboard 

1 

Vertical  Bending 

4.31 

4.31 

4.31 

Nacelle 

Coupled  Side 

Bending  &  Torsion 

2.07 

2.07 

2.07 

Outboard 

1 

Vertical  Bending 

4.28 

4.28 

4.28 

Nacelle 

Coupled  Side 

Bending  &  Torsion 

2.24 

2.24 

2.24 

*Includes  horizontal  and  vertical  tail  masses. 

**See  Fig.  23 


TABLE  7 


REPRESENTATIVE  FREE-FREE  NORMAL  MODES 
VERTICAL  ANALYSES 


Mode 

No. 

Responding 

Component 

Modal  Frequency, 

cps 

Structural 

Damping, 

g** 

OWE 

Half 

Fuel 

Full 

Fuel 

Ground 

Vibration 

Test* 

1 

Wing  Bending 

1.40 

1.35 

1.08 

1.44 

•  040 

2 

Nacelle 

2.07 

2.07 

2.05 

2.07 

.027 

3 

Nacelle 

2.22 

2.20 

2.14 

— 

.037 

i* 

Wing  Torsion 

2.46 

2.44 

2.25 

2.96 

.022 

5 

Forward  Fuselage 

3.32 

3.26 

3.19 

3-75 

.030 

6 

Fuselage 

4.07 

3-65 

3.40 

5.40 

.030 

7 

Wing  Bending 

6.07 

5.30 

4.55 

6.50 

.030 

*The  ground  vibration  teat  was  conducted  for  the  720  airplane  with 
the  main  gear  and  nose  gear  restrained.  (Ref.  17) 


**Damping  coefficient,  g,  is  equal  to  twice  the  actual  damping 
divided  by  the  critical  damping. 


TABLE  8 

CANTILEVERED  MODES  LATERAL  ANALYSES 


Airplane 

Component 

No.  of 
Discrete 
Masses* 

Cantilevered 

Mode 

Description 

Modal 

Frequency, 

cps 

Aft 

Fuselage 

12 

1st  Side 

Bending 

2.U1 

Aft 

Fuselage 

12 

2nd  Side 

Bending 

7.22 

Aft 

Fuselage 

12 

1st  Torsion 

5.83 

Forward 

Fuselage 

7 

1st  Side 

Bending 

3.66 

Fin 

13 

1st  Bending 

4.3!* 

Fin 

13 

1st  Torsion 

14.62 

**T<*m- 


Equations  of  Motion 

The  equations  of  motion  for  the  720B  gust  analyses  were  developed  by 
Lagrange's  energy  approach.  The  derivation  of  vertical  and  lateral 
equations  of  motion  is  discussed  in  Appendix  A.  The  equations  for  both 
analyses  are  a  set  of  simultaneous  second-order  differential  equation 
in  the  generalized  coordinates  described  previously.  These  equations 
can  be  expressed  in  matrix  form  as  follows : 


D.  . 
id 
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a. 
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C .  a 
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The  coefficient  matrices  of  q.^,  q^,  and  qi  are  the  generalized  mass, 

.[Mij]  >  the  structural  and  aerodynamic  damping,  [D^]  ,  and  the  structural 
and  aerodynamic  stiffness,  [K^j]  ,  respectively.  The  right  hand  side  of 
the  equation  represents  the  generalized  forces  in  the  degrees  of  freedom 
due  to  the  gust  excitation,  cJ  . 

The  generation  of  the  generalized  mass  coefficient  matrix  began  with  the 
lumped  mass  representation  of  the  airplane  in  terms  of  discrete  masses 
and  mass  moments  of  inertia  about  the  discrete  mass  center  of  gravity 
for  each  mass  panel.  Then,  the  center  of  gravity  mass  properties  were 
transferred  to  their  respective  elastic  axis  reference  stations.  The 
total  kinetic  energy  of  the  system  was  written  in  terms  of  generalized 
coordinate  velocities.  Differentiation  of  the  kinetic  energy  expression, 
first  with  respect  to  the  generalized  coordinate  velocities,  and  then 
with  respect  to  time,  resulted  in  the  generalized  mass  matrix, 

The  damping  coefficients  arise  partly  from  the  aerodynamic  response 
forces  and  partly  from  structural  damping.  The  aerodynamic  damping  re¬ 
sults  from  that  portion  of  the  total  aerodynamic  force  which  is  in  phase 
with  the  particular  generalized  coordinate  velocities.  The  structural 
damping  is  represented  as  forces  that  are  proportional  to  the  structural 
stiffnesses  but  are  in  phase  with  the  response  velocities.  The  values 
for  sti'uctural  damping  used  in  the  vertical  analyses  were  obtained  from 
ground  vibration  tests  of  the  r'20  airplane.  The  damping  coefficient  for 
each  mode  for  the  vertical  analyses  is  listed  in  Table  7*  The  energy  or 
work  expression  for  these  damping  forces  was  differentiated  with  re¬ 
spect  to  each  coordinate  to  arrive  at  the  generalized  damping  matrix, 

[»«]• 

The  stiffness  matrix  also  results  from  structural  and  aerodynamic 
sources.  The  stiffness  of  the  primary  airplane  structure  is  represented 
in  the  equations  of  motion  by  the  elastic  potential  energy  stored  in  the 
structure  when  it  is  deformed  in  a  linear  combination  of  the  selected 
vibration  mode  shapes.  The  aerodynamic  work  resulting  from  the  response 


displacements  of  the  primary  components  of  the  airplane  is  added  to  the 
internal  elastic  energy  and  the  entire  expression  is  differentiated  with 
respect  to  each  generalized  coordinate  to  form  the  stiffness  matrix, 

[Kij]  ' 

The  yaw  damper  effects  in  the  lateral  equations  of  motion  are  handled  by 
including  a  generalized  force  specifically  for  the  yaw  damper.  This 
generalized  force  is  a  function  of  both  the  gust  excitation  frequency 
and  the  resulting  airplane  yawing  velocity.  The  yaw  damper  system  con¬ 
sists  of  a  yaw  rate  sensor,  a  yaw  axis  amplifier,  a  servo  system  with 
rudder  position  and  rudder  angular  rate  feedback  loops,  and  an  hydraulic 
rudder  actuator  unit. 

In  the  airplane,  the  signal  from  the  yaw  rate  sensor  is  filtered  to  re¬ 
move  the  effects  of  both  slow  maneuvers  and  higher  frequency  fuselage 
vibrations.  The  filtered  signal  is  amplified  and  applied  to  an  electronic 
servo  which,  in  turn,  drives  the  hydraulic  rudder  actuator. 

The  yaw  damper  system  is  simulated  in  the  lateral  analyses  by  incorpora¬ 
ting  a  transfer  function  due  to  yawing  velocity.  This  transfer  function 
is  a  composite  of  the  transfer  functions  for  all  of  the  sub-parts  of  the 
yaw  damper  system.  The  generalized  force  associated  with  the  rudder 
deflection  varies  with  frequency,  but  it  is  a  function  of  yaw  rate  only 
at  any  given  frequency.  The  generalized  force  coefficients  were  added 
to  the  coefficients  in  the  rudder-fixed  equations  of  motion  to  account 
for  the  yaw  damper  effects.  A  more  complete  description  of  the  mathe¬ 
matical  model  for  the  ya w  damper  is  given  in  Appendix  A,  Equations  of 
Motion. 

Both  the  vertical  and  lateral  equations  of  motion  were  solved  to  obtain 
the  real  and  imaginary  parts  of  the  steady-state  complex  responses  of  the 
generalized  coordinates  due  to  a  1.0  ft.  per  sec.  continuous  sinusoidal 
gust  velocity  excitation.  Consequently,  the  accelerations,  q^,  could 
be  replaced  by  ~tu£q^  and  the  velocity,  q^,  by  iruq^,  where  <y  is  the  ex¬ 
citation  frequency.  When  these  assumptions  were  incorporated  into  the 
equations  of  motion,  the  matrix  coefficients  ,  ic^D-^]  and 

could  be  combined  into  a  single  complex  set  of  matrix  coefficients 
[Aij  +  iB^ J  of  the  generalized  coordinates  for  each  frequency.  The 
generalized  coordinate  responses  were  obtained  by  premultiplying  the 
generalized  gust  forces,  j  CM  ,  by  the  inverse  of  this  matrix  for  each  of 
the  selected  frequencies  '  ' 

Solutions  were  obtained  in  the  vertical  analyses  to  define  the  complex 
frequency  response  functions  over  the  range  from  0. 5(lO)“5  to  0.6  rad. 
per  ft.  The  sensitivity  of  the  rms  values,  A,  and  the  number  of  zero 
crossings,  N0,  to  the  upper  and  lower  limits  of  the  frequency  range  and 
to  the  overall  number  of  frequencies  was  examined.  The  lower  end  of  the 
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frequency  range  was  taken  at  0.  5(lO)“5  rad.  per  ft  which  was  well  below 
the  short  period  mode  response  frequency.  The  maximum  frequency,  0.6  rad. 
per  ft. ,  was  sufficiently  high  such  that  the  number  of  zero  crossings 
had  converged.  A  total  of  95  frequencies  between  the  above  limits  were 
used  in  the  vertical  analysis.  The  major  contribution  to  the  rms  values 
is  in  the  lower  end  of  the  frequency  range;  therefore,  the  solution  fre¬ 
quencies  were  selected  at  closer  intervals  in  the  lower  frequency  range. 

Load  Equations 


The  detailed  loads  for  the  aircraft  structure  were  obtained  from  the 
complex  frequency  responses  of  the  generalized  coordinates.  This  was 
accomplished  by  transforming  the  generalized  inertia,  damping  and 
stiffness  forces  uscu  in  the  equations  of  motion  into  shears,  bending 
moments,  and  torsions  referenced  to  given  stations  on  the  elastic  axes. 
The  shear,  moment,  and  torsio*-  coefficients  were  calculated  for  unit 
deflections  in  the  generalizeu  coordinates.  Then,  these  coefficients 
were  multiplied  by  the  complex  frequency  responses  of  the  generalized 
coordinates  to  obtain  load  complex  frequency  responses  for  all  elastic 
axis  reference  stations. 

Complex  frequency  response  functions  for  stresses  were  obtained  by  in¬ 
corporating  unit  load  solution  stress  coefficients  from  the  airplane 
stress  analysesv11^l2A3)  into  the  load  equations. 

Spectral  Analysis 


The  complex  frequency  response  functions  for  the  incremental  shears, 
moments,  and  torsions  were  used  to  define  complex  frequency  response 
functions  for  incremental  axial  and  shear  stresses  in  the  detailed  air¬ 
plane  structure.  The  output  power  spectra  for  the  loads  and  stresses 
were  obtained  as  follows: 


*0(n>  =*v(n) 


T(iO) 


(5) 


where, 


output,  load  or  stress  power  spectral  density  function 

(ft)  input,  normalized  vertical  or  lateral  gust  power  spec¬ 
tral  density  function 

T(ift)  complex  frequency  response  function  for  load  or  stress 
at  a  given  station  in  the  airplane 


40 


The  rms  values  for  the  loads  and  stresses  for  a  rms  gust  velocity  of 
unity,  A,  were  determined  by  integrating  the  appropriate  output  power 
spectrum. 
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The  zero  crossings  per  ft.  for  all  loads  and  stresses  were  determined 
by  Rice's  expression'^'. 
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Since  the  actual  strength  of  a  detailed  structural  element  is  in  general 
governed  by  the  combined  stresses,  stress  interaction  functions  and  one 
factor  level  flight  stresses  must  be  incorporated  into  the  statistical 
analysis.  The  problem  was  to  determine  the  probability  of  exceeding  the 
limit  strength  for  a  given  flight  condition  (Design  Envelope  Approach) 
and  to  determine  the  number  of  times  the  limit  strength  is  exceeded  per 
unit  time  (Flight  Profile  Approach).  The  joint  probability  density 
function  for  two  components  of  stress  on  a  structural  element,  say,  axial 
stress  and  shear  stress  is, 
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where, 


axial  stress,  psi 
shear  stress,  psi 


one  factor  axial  and  shear  stresses,  respectively,  psi 
rms  axial  and  shear  stresses,  respectively,  psi 
correlation  coefficient 
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Statistical  Analysis  -  Design  Envelope  Approach 


The  volume  under  the  probability  density  functionshown  in  Fig.  14  bounded 
by  the  interaction  limit  strength  enve-lope  is  equal  to  the  probability 
that  the  limit  strength  will  not  be  exceeded,  that  is,  that  the  margin 
of  safety  is  greater  than  zero.  The  probability  that  the  limit  strength 
will  be  exceeded  is  1.0  minus  the  probability  that  the  limit  strength 
will  not  be  exceeded.  Therefore, 


P(MS  <  0)  =  1  -  P(MS  ^  0) 


(9) 


Since  the  input  or  gust  power2spectral  density  varies  directly  with  the 
mean  squared  gust  velocity,  or  Oy  ,  the  probability  of  exceeding  the 
limit  design  strength  is  a  function  of  aw  or  That  is, 
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Consider  a  stiffened  skin-stringer  segment  typical  of  those  used  in  wing 
panel  construction  of  subsonic  jet  transport  airplanes.  A  cross  section 
through  the  segment  normal  to  the  stringer  is  shown  in  Fig.  15-  It  was 
assumed  that  the  airplane,  of  which  this  structural  element  is  a  part, 
was  flying  in  turbulent  air.  It  was  of  interest  to  determine  (l)  the 
probability  of  exceeding  the  limit  design  strength  of  the  segment  and 
(2)  the  expected  number  of  times  per  hour  that  the  limit  strength  would 
be  exceeded  for  various  root  mean  square  (rms)  gust  velocities,  o^< 

Numerical  information  pertinent  to  this  example  is  given  in  Table  9*  The 
appropriate  limit  strength  interaction  diagram  is  shown  with  a  sketch  of 
the  structural  segment  in  Fig.  15*  It  will  be  noted  that  two  level- 
flight  stress  points,  skin  stress  and  segment  stress,  are  shown  in 


TABLE  9 

NUMERICAL  EXAMPLE 
STRESS  DATA  FOR  STATISTICAL  ANALYSIS* 

Level  Flight  Stresses  (Based  on  Appropriate  Effective  Areas) 

Skin  Stress,  Tension  f0t,  psi  15,600 

Segment  Stress,  Compression,  f0  ,  psi  15,200 

Shear,  £0,  psi  c  -2,190 

Dynamic  Analysis  Stress  Data  (Power  Spectral  Analysis- RMS  Gust  Velocity, 
0  =  1.0  ft.  per  sec.  ) 

RMS  Tension  Stress,  Af^.,_psi  184.19 

RMS  Compressive  Stress,  Ap  ,  psi  179- 51 

RMS  Shear  Stress,  A^  ,  psiC  38.53 

Correlation  Coefficient  -0.2141 

(Axial  and  Shear  Stresses) 

Zero  Crossings  per  sec 

Axial  Stress,  N0„  1.2205 

Shear  Stress,  N0£  2.7235 

RMS  Stress  Rates,  psi  per  ft .per  sec 

Tension,  CT^  1,413 

Compression,  crpc  1,377 

Shear,  cra  659 

Allowable  Limit  Design  Stresses,  psi 

Tension,  F^  40, 000 

Compression,  Fc  21,270 

Shear,  Ft  20, 500 


*See  Appendix  B  for  development  of  methods  for  Statistical  Analysis 
of  Combined  Random  Stresses 


Table  9  and  in  Fig.  15*  These  are  reference  points  for  two  joint  pro¬ 
bability  stress  distributions.  The  compression-shear  or  segment  stress 
probability  distribution  was  applied  to  the  compression  portion  of  the 
interaction  diagram  in  Fig.  15,  and  the  tension-shear  or  skin  stress 
distribution  was  applied  to  the  tension  portion  of  the  diagram.  The 
volume  of  the  compression-shear  joint  probability  density  function 
within  the  compression  region  of  the  interaction  diagram  was  added  to 
the  volume  of  the  tension-shear  density  function  in  the  tension  region 
to  determine  the  probability  that  the  limit  design  strength  will  not  be 
exceeded.  A  typical  plot  of  ?(MS  <  0 ,aw)  versus  <7W  is  shown  in  Fig.  l6. 

Certain  structural  elements  in  the  fuselage  near  either  the  vertical  or 
lateral  bending  axis  are  only  subjected  to  fluctuating  compression 
stresses  and  shear.  They  cannot  be  stressed  in  tension  as  a  result  of 
the  gust  loading.  This  occurs,  because  the  location  of  the  neutral 
plane  through  the  fuselage  section  shifts  across  the  element  as  the 
direction  of  the  bendipg  moment  on  the  fuselage  changes  sign.  The 
element  is  always  on  the  compression  side  of  the  neutral  plane.  This 
compression-compression  situation  was  handled  in  the  statistical  analy¬ 
sis  of  limit  strength  exceedances  by  removing  the  tension  part  of  the 
interaction  diagram  and  by  substituting  another  compression,  part.  Now, 
these  special  cases  could  be  treated  in  the  same  manner  as  normal 
elements. 

The  objective  of  the  design  envelope  approach  was  to  determine  a  level 
of  gust  i  tensity,  which  when  multiplied  by  a  stress  A  value  for 

a  critical  element  would  just  equal  the  limit  design  strength  for  the 
element.  Thus,  the  design  envelope  approach  implies  a  one-dimensional 
probability  distribution  for  design  stress.  The  statistical  distribu¬ 
tion  of  stress  values  is  assumed  to  be  Gaussian  and  can  be  completely 
described  by  the  rms  stress  value  for  a  one-dimensional  distribution. 

When  the  probability  of  exceeding  limit  strength  is  governed  by  a  joint 
probability  function  in  two  dimensions,  three  parameters,  i.e.  two  rms 
stress  values  and  a  correlation  coefficient,  are  required.  As  a  result, 
no  single  statistical  parameter  can  be  scaled-up  to  a  design  strength 
level.  However,  the  probability  of  exceeding  the  design  strength  of  the 
element  can  be  matched  to  i  probability  for  exceeding  an  incremental 
eg  acceleration*  for  the  same  flight  condition.  The  distribution  of 


The  development  that  follows  is  essentially  unchanged  if  the  c.g.  ac¬ 
celeration  is  replaced  by  a  generalized  load  quantity,  y/A.  Thus  it 
can  be  seen  that  the  numerical  results  are  not  dependent  upon  the  choice 
of  the  particular  output  quantity  for  which  the  single-parameter  proba¬ 
bility  distribution  is  obtained. 
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incremental  eg  acceleration  is  a  single  parameter  distribution, 
unique  can  be  determined  for  that  level  of  incremental  eg 
acceleration. 

The  probability  of  exceeding  specified  values  of  incremental  eg  accel¬ 
eration  is  determined  as  follows: 
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This  equation  is  evaluated  for  various  specified  A a  values  to  obtain 
the  distribution  functions  shown  in  Fig.  17»  It  will  be  noted  that  the 
diagrams  in  Fig.  17  are  universal,  and  can  be  applied  to  any  airplane 
at  any  altitude. 

Of  course,  the  cumulative  distribution  functions  shown  in  Figs.  1 6  and  17 
will  diff  r,  because  the  first  is  based  on  joint  probability  considera¬ 
tions,  an-  the  second  is  a  single  parameter  distribution.  The  two  dis- 
tributi  ..  "'unctions  cannot  be  made  to  coincide  everywhere.  The  question 
was,  where  s.  nxld  they  be  matched? 

When  each  probability  function  is  multiplied  by  f(crw),  the  probability 
density  of  rms  gust  velocity,  and  integrated  over  all  ow,  we  obtain  a 
new  probability  number  indicating  the  chance  of  exceeding  the  design 
strength  at  any  instant  in  time  for  all  weather  and  seasonal  conditions. 
Typical  overall  probability  density  functions  for  incremental  eg  accel¬ 
eration  are  shown  in  Fig.  18.  Similar  typical  information  for  a  struc¬ 
tural  element  is  shown  in  Fig.  l$h 

It  is  ovbious  that  one  curve  of  the  family  of  curves  shown  in  Fig.  18 
will  best  match  the  curve  in  Fig.  19-  The  eg  acceleration  curve  that 
best  matches  the  curve  for  the  structural  element  is  the  one  having  very 
nearly  the  same  area  when  integrated.  It  is  noted  that  the  maxima  of 
the  overall  probability  density  curves  for  eg  acceleration  lie  on  a 
straight  line  as  shown  in  Fig,  18,  and  it  is  possible  to  cross  plot  the 
maximum  overall  probability  density  versus  acceleration  level  as  shown 
in  Fig.  20.  This  diagram  is  unique  for  the  flightAcondition.  The 
straight  line  will  intersect  the  vertical  axis  at  f (aw)  =  r(0)  since 
P(An  >  0)  =  1.0.  Now  it  is  possible  to  determine  for  the  struc¬ 
tural  element  by  determining  the  maximum  probability  density  in  Fig.  19, 
(9. 0(l0) _13).  This  maximum  probability  density  was  used  in  Fig.  20  to 
determine  An^,  (1.89),  and  for  the  element. 
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FIGURE  18.  OVERALL  PROBABILITY  DENSITIES  FOR  EXCEEDING 
SPECIFIC  INCREMENTAL  CG  ACCELERATIONS 
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It  is  obvious  that  the  limit  design  gust  intensity,  o^rj^}  would  vary 
with  altitude,  since  the  statistical  characteristics  of  atmospheric 
turbulence  changes  with  altitude.  This  variation  is  accounted  for  by 
use  of  Eq.  (l4). 

-y/b,A  -y/bpA 

N(y)=F1N0e  +PgN0e  (l4) 


When  we  consider  y/A  to  be  a  limit  design  gust  velocity,  <7W%>  then, 
Eq.  (l4)  is  rewritten  as  follows: 


N(<V?d)  -K%)/bi 

— ?? - =P.  e 


-(<V7a)/b2 

+  P2  e 


(A  =  1.0) 


(15) 


The  value,  N(crw7/^)/N0,  is  thought  of  as  an  exceedance  ratio.  It  is  the 
ratio  of  the  number  of  gusts  exceeding  the  limit  design  gust  per  unit 
time  or  distance  at  a  given  altitude  to  the  total  number  of  gusts  per 
unit  time  or  distance  at  that  altitude. 

The  analysis  procedure  for  the  Design  Envelope  Approach  was  as  follows: 
(l)  Determine  the  lowest  allowable  values  of  <7wJfa  or  <7V for  the  major 
components  of  the  reference  airplane,  the  720B,  (2)  Determine  the  ex¬ 
ceedance  ratio  for  the  critical  elements  or  components  of  the  airplane 
for  a  particular  critical  altitude,  (3)  Establish  an  acceptable  ex¬ 
ceedance  ratio  for  each  altitude  from  all  design  envelope  analysis 
conditions.  Determine  the  limit  design  aw%  for  each  altitude. 

Statistical  Analysis  -  Flight  Prefile  Approach 


The  Flight  Profile  Approach  involves  the  determination  of  the  average 
number  of  times  per  hour  that  the  limit  strengths  of  the  mere  critical 
structural  components  or  elements  are  exceeded.  The  average  number  of 
limit  strength  exceedances  is  a  weighted  average  number  derived  from  the 
estimated  airline  usage  of  the  (20B  airplane. 
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The  average  number  of  limit  strength  exceedances  per  hour  for  a  given 
flight  condition  is  established  by  determining  the  number  ;of  times  per 
foot  travelled  that  the  limit  strength  envelope-  is  perforated  or  crossed, 
Nc«-  It  -is  assumed  that  the  number  of  limit  strength  exceedances  is 
equal  to  one-half  the  number  of  perforations  of  the  strength  envelope. 
That  is,  there  is  one  peak  value  exceeding  limit  strength  for  each  two 
crossings.  One-half  of  the  number  of  strength-  envelope  crossings  per 
foot  travelled  times  -the  airplane  velocity  in  ft.  per  sec.  times  3,600 
sec.  per  hour  gives  the  number  of  limit  strength  exceedances  per  flight 
hour  for  a  given  rms  gust  intensity  oy  or  <7V.  A  typical  diagram  show¬ 
ing  the  limit  strength  exceedances  per  ft.  travelled  for  various  rms 
gust  intensities  is  shown  in  Fig.  21. 

In  order  to.-  account  for  the  variation  in  rms  gust  intensity,  the  limit 
strength  exceedances  per  hour  is  multiplied  by  the  probability  density 
function  for  rms  gust  velocity  as  shown  in  Figj_  22.  The  integral  of 
this  new  probability  function  over  all  (rw  is,  G,.  the  expected  number  of 
limit  strength  exceedances  per  hour.  That  is, 
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The  limit  strength  exceedances  per  hour  for  the  representative  flight 
profiles  is  obtained  by  multiplying, the  exceedances  per  hour  for  each 
flight  condition  in  the  profile  by  the  fraction  of  the  total  profile 
time  spent  at  that  condition,  and  then  summing  the  exceedances  per  hour 
for  all  conditions  in  the  particular  profile. 

The  'exceedances  per  average  flight  hour  for  'total  airplane  usage  is 
then-  determined  by  weighting  the  exceedances  per  hour  of  each  profile 
by  the  fraction  of  total  airplane  usage  represented  by  the  particular 
profile:  and  then  summing  over  all  the  profiles.  This  process  is  re¬ 
presented  by  the  following  equation: 
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where 

Np  total  airplane  exceedances  per  average  flight  hour 

G.  exceedances  per  hour  for  a  particular  flight  condition 
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Nc/2,  EXCEEDANCES  OF  LIMIT  DESIGN  STRENGTH  PER  FT 


EXCEEDANCES  OF  LIMIT  DESIGN  STRENGTH 


DISTRIBUTIONS  OF  LIMIT  DESIGN  STRENGTH 


<rw,  RMS  GUST  VELOCITY  (ft/s*c) 

FIGURE  22.  DISTRIBUTIONS  OF  LIMIT  DESIGN  STRENGTH 
EXCEEDANCES  WITH  RMS  GUST  INTENSITY 
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the  time  spent  at  a  particular  .flight  condition 
the  total  profile  flight  time 

fraction  of  total  airplane  usage  represented  by  a 
particular  profile 

the  number  of  conditions  in  a  particular  profile 

the  number  of  profiles  used, to  represent  the  airplane 
usage 


RESULTS 


General  -Remarks 


In  the  past,  airplane  .power  spectral  gust  analyses  results  have  been 
used  to  determine  the  effects  of '  airplane  dynamic  responses  on  the 
bending  moments,  .shears,  and  torsions-  for  the  major  components  of  an 
airplane;  These  incremental  loads,  usually  the  wing  bending  moments , 
were  compared  with  comparable,  loads  from  a  static  aeroelastic  analysis, 
using  the  gust  load  formula,  to  arrive  at  a  dynamic  factor.  Then,  the 
■dynamic  factor  was-  applied  to  other  aeroelastic  solutions  to  obtain  de¬ 
sign  incremental  loads  for  all  gust,  design  conditions. 

The  basic  results  of  a  power  spectral  analysis,  the  rms  loads,  A,  and  the 
number  of  times  these  loads  cross  their  mean  values  with  positive-  (or 
negative)  slope  per  sec.  or  per  ft;,  N0,  are  not  in  themselves,  adequate 
to  assess  airplane  strength.  However,  the  bending- moments,  shears,  and 
torsions,-  and  the  incremental  eg  accelerations  do  provide  insight  into 
an  airplane  dynamic  responses  in  rough  air,  and  can  be  used  as  guides  to 
ascertain  critical  locations  in  the  airplane  structure.  Then,  the  prop¬ 
erly  phased  combined  effects  of  the  bending  moments,  shears,  and  torsions 
should  be  used-  to  assess  the  strength  at  these  critical  locations. 

The  rms  bending  moments,  shears,  and  torsions  for  the  720B  airplane  were 
determined  at  the  locations  indicated  in  Figs.  23  and  24  for  the  vertical 
and  lateral  analyses,  respectively.  The  combined  effects  of  vertical 
and  lateral  gust  loads  were  studied  for  two  locations  on  the  aft  fuse¬ 
lage,  Body  Balance  Stations  1000  and  1280.  Stress  analyses  were  con¬ 
ducted  for  the  most  critical  locations  on  the  wing,  fuselage,  and  the 
vertical  tail.  These  locations  consisted  of  one  wing  location,,  Wing 
Station  360  at  Body  Buttock  Line  260  or  33  percent  of  the  wing  semi-span 
(  1-  *33) •  The  stresses,  were  studied  for  three  body  locations,  Body 

Balance  Stations  480,  1000,  and  1280.  One  location  on  the  vertical  tail, 
Fin  Station  1 65.5,  Elastic  Axis  (EA)  .Station  158  was  investigated. 


Correlation  of  Dynamic  and  Aeroelastic  Analyses 

The  total  loads  and  stresses  that  were  used  to  assess  the  airplane 
strength  resulted  from  the  incremental  dynamic  loads  and  the  one-factor 
level  flight  loads.  The  level  flight  loads  were  obtained  from  aero¬ 
elastic  solutions, (6*7)  and  the  incremental  dynamic  loads  were  obtained 
from  power  spectral  analysis.  In  order  to  determine  if  these  two  inde¬ 
pendent  solutions  represented  the  airplane  in  a  consistent  manner,  the 
one-factor  level  flight  loads  were  obtained  for  one  analysis  condition, 
Condition  1,  using  the  dynamic  analysis  equations  of  motion.  This  was 
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accomplished  by  setting  the  pitch  and  flexible  generalized  coordinate 
accelerations,,  the  pitch  and  translation  displacements,  and  all  of  the 
generalized  coordinate  velocities  equal  to  zero.  The  vertical  trans¬ 
lation  acceleration  was  set  at  l.Og.  In  order  to  allow  for  a  force 
balance  in  the  vertical  translation  equation  and  a  moment  balance  in  the 
pitch  equation,  two  unknowns  were  introduced  to  replace  the  airplane 
pitch  and  translation  displacement  freedoms.  One  of  the  Unknowns  was 
the  -angle  of  attack  which  gives  rise  to  the  lift  necessary  to  offset  the 
airplane  weight  plus  the  balancing  tail  load.  The  other  unknown  was  the 
horizontal  tail  load  which' must  balance  the  airplane  pitching  moments. 

The  wing  bending  moments,  shears,  and  torsions*  and  the  wing  static 
deflections  resulting  from  the  reduced  dynamic  solution  were  obtained  to 
compare  with  corresponding  values  from  the  aeroelastic,  solution.  The 
aeroelastic  wing  loads  shown  in  this  comparison  were  obtained  by  sub¬ 
tracting  a  zero  load  factor  solution  from  a  one-factor  solution.  The 
comparison  is  made  on  the  basis  of  a  static  incremental  one-factor  load. 
The  results  in  Fig.  25  show  the  wing  bending  deflections  obtained  from 
the  two  solutions.  The  wing  bending  moment,  shear,  and  torsion  loads 
are  compared  in  Fig.  :2 6.  The  single  of  attack  .required  for  level  flight 
was  4.27  degrees  for  the  reduced  dynamic  solution  and  4.32  for  the  aero¬ 
elastic  solution. 


Frequency  Response  Functions 

The>  .equations  of  motion  were  solved  tp  obtain  steady  state  complex  fre¬ 
quency  response  functions  for  the  generalized  coordinates  due  to  1.0  ft. 
per  sec.  sinusoidal  gust  excitation.  Examples  of  the  modulus  of  these 
frequency  response  functions  for  the  Vertical  Analyses,  Condition  24, 
are  presented  in  Figs.  27(a)  through  27(j).  The  response  quantities  are 
identified  as  qj_  through  q^  for  the  seven  free-free  airplane  symmetrical 
modes  and  q^,  qp  and  qx  for  the  rigid  airplane  freedoms  of  vertical  trans¬ 
lation,  pitch,  and  fore  and  aft  translation,  respectively.  These  fre¬ 
quency  response  functions  illustrate  the  predominate  responses  for  each 
individual  generalized  coordinate,  and  indicate  the  coupling  that  exists 
between  the  coordinates.  The  frequency  response  curve  for  airplane 
pitch,  Fig.  27(i),  shows  the  well  damped  airplane  short  period  or  pitch 
mode  at  O.36  cps;  0.028  rad.  per  ft.  The  airplane  pitch  response  has  a 
large  effect  on  the  airplane  loads,  because  of  the  high  energy  in 
atmospheric  turbulence  within  this  frequency  range. 

Typical  generalized  coordinate  response  curves  for  the  lateral  analyses, 
Conditions  10  and  10YD,  are  shown  in  Fig.  28(a)  through  28(i)  both  with 
and  without  the  yaw  damper  operating.  The  Dutch  roll  response  at  0.20 
cps  or  0.0015  rad.  per  ft.  is  the  major  contributor  to  the  lateral  loads. 
Thus,  the  changes  which  effect  the  Dutch  roll  response  have  a  large  effect 
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27a.  FREQUENCY  RESPONSE,  FIRST  SYMMETRICAL 
AIRPLANE  MODE 


FIGURE  27b.  FREQUENCY  RESPONSE,  SECOND  SYMMETRICAL 
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FIGURE  27c.  FREQUENCY  RESPONSE,  THIRD  SYMMETRICAL 

AIRPLANE  MODE 
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FIGURE  27 d.  FREQUENCY  RESPONSE,  FOURTH  SYMMETRICAL 

AIRPLANE  MODE 


FIGURE  27e.  FREQUENCY  RESPONSE,  FIFTH  SYMMETRICAL 

AIRPLANE  MODE 
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FIGURE  27 f.  FREQUENCY  RESPONSE,  SIXTH  SYMMETRICAL 

AIRPLANE  MODE 
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FIGURE  28d.  FREQUENCY 
RESPONSE,  FIRST 
CANTILEVERED  AFT  BODY 
SIDE  BENDING  MODE 


NOTE: 


YAW  DAMPER  ON 
YAW  DAMPER  OFF 
CONDITION  10,  10 


FIGURE  28e.  FREQUENCY  RESPONSE, 
SECOND  CANTILEVERED  AFT  BODY 
SIDE  BENDING  MODE 


YD 


0  (rad/Tt) 


FIGURE  28f.  FREQUENCY  RESPONSE, 
FIRST  CANTILEVERED  AFT  BODY 
TORSION  MODE 
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FIGURE  28g.  FREQUENCY 
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FIGURE  28h.  FREQUENCY 
RESPONSE,  FIRST 


FIGURE  28i.  FREQUENCY 
RESPONSE,  FIRST 
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on  magnitudes  of  the  lateral  loads.  When  the  yaw  damper  is  engaged,  the 
rudder  is  commanded  by  the  yaw  damper  servo  system  and  will  deflect  in 
a  manner  to  provide  positive  damping  in  the  Dutch  roll  mode.  The  rer 
suits  of  the  yaw  damper  operation  is  readily  apparent.  By  observing 
the  decrease  in  the  rigid'  airplane  yaw  response  not  only  are  the  fre¬ 
quency  response  functions  for  the  rigid,  airplane  motions  attenuated,  but 
the  frequency  at  .which  the  maximum  response  occurs  is  reduced.  This  is 
reasonable,  because  at  the  Dutch  roll  frequency,  0.07  cps,  the  phase 
shift  through  the  yaw  damper  system  causes  the  rudder  displacement  to 
lead  the  yaw  rate  signal.  This  results  in  a  rudder  deflection  that  re¬ 
duces  the  fin  aerodynamic  force  due  to  yawing  displacement  which,  in 
turn,  causes  a  reduced  Dutch  roll  resonant  frequency. 

The  coupling  between  the  flexible  modes  and  the  Dutch  roll  mode  occurs 
■at  a,  higher  frequency  with  the  autopilot  engaged.  This  may  seem  to 
contradict  the  effect  evidenced  by  the  rigid  airplane  response;  however, 
this  difference  does  exist,  and  there  is  a  logical  explanation.  The 
explanation  again  lies  with  the  phasing  of  the  rudder  deflection  caused 
by  the  yaw  damper  output.  In  this  case;  the  frequency,  .63  cps,  0.0046 
rad.  per  ft.,,  is  such  that  the  rudder  displacement  lags  the  yaw  rate 
signal;  therefore,  the  rudder  deflection  adds,  to  the  fin  aerodynamic 
force.  This  increase  in  force  results  in  a  larger  flexible  mode  excita¬ 
tion  at  this  frequency.  As  the  response  frequency  increases  toward  a 
structural  resonance  frequency  the  cut-off  filter  in  the  yaw  damper 
amplifier  system  attenuates  the  yaw  damper  output  so  that  self-excitation 
will  not  occur. 

Complex  frequency  responses  for  the  loads  were  obtained'  by  substituting 
the  complex  generalized  coordinate  responses;  into  the  airplane  load 
equations.  These  load  equations  sum  the  incremental  inertia  loads,  the 
airloads  resulting  from  the  coordinate  responses,  and  the  applied  gust 
loading.  Fig.  29  shows  the  load  responses  for  Vertical  Analyses,  Condi¬ 
tion  24,  for  wing  bending  moment,  shear,  and  torsion  at  the  most  critical 
wing  station,  Wing  Eta  Station  0.33*  Similar  load  frequency  response 
functions  are  shown  for  Wing  Eta  Station  0.12.  Fuselage  vertical  shear 
and  bending  moment  frequency  responses  at  the  critical  forward  body  sta¬ 
tion,  Body  Station  540  and  at  Aft  Body  Station  1360,  are  shown  in 
Fig.  30. 

Fig.  31  gives  the  vertical  tail  load  frequency  responses  for  lateral 
analysis  Conditions  10  and  10YD.  These  frequency  responses  are  for 
vertical  tail  bending  moment,  shear,  and  torsion  at  the  vertical  tail 
Elastic  Axis  Stations  67  and  158  with  and  without  the  yaw  damper  oper¬ 
ating,  Fig.  32  shows  the  lateral  bending  moment,  shear,  arid  torsion 
frequency  response  functions  for  Body  Station  1375  with  and  without  the 
yaw  damper  operating. 
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FIGURE  30a.  FREQUENCY  RESPONSE, 
FUSELAGE  VERTICAL  SHEAR 
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FIGURE  30b.  FREQUENCY  RESPONSE,  FUSELAGE  VERTICAL 
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FIGURE  30d,  FREQUENCY  RESPONSE,  FUSELAGE  VERTICAL 
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FUSELAGE  SHEAR  (1<T  Ib/Ff/tec) 


FIGURE  32a.  FREQUENCY 
RESPONSE,  FUSELAGE 
LATERAL  SHEAR 


NOTE 
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CONDITION  10,  10  YD 
BODY  STATION  1375 


FIGURE  32b.  FREQUENCY 
RESPONSE,  FUSELAGE 
LATERAL  BENDING  MOMENT 


RESPONSE,  FUSELAGE 
TORSIONAL  MOMENT 


The  major  response  peaks  for  the  wing  loads  result  from  airplane  pitch 
and  from  the  first  and  fifth  flexible  modes.  This  is  reasonable,  be¬ 
cause  the  first  and  fifth  flexible  modes  constitute  the  largest  part  of  j 

the  coupled  wing. bending  deflection.  The  wing  torsion  loads  receive  ' 

their  major  responses  from  the  airplane  pitch  and  first  and  fourth 
flexible  modes,  the  fourth  mode  is  primarily  wing  torsion.  The  second  r 

and  third  flexible  modes  are  predominantly  nacelle  side  bending  modes;  [ 

therefore,  they  contribute  relatively  little,  to  the  overall  wing  loads. 

The  fuselage  vertical  bending  moment  and  shear  loads  are  derived  mainly 

from  the  fifth  flexible  airplane  mode,  which  is  predominantly  fuselage  If 

vertical  bending.  * 


The  primary  response  for  the  vertical  tail  results  from  the  Dutch  roll 
mode;  however,,  when  the  yaw  damper  is  on  the  fin  and  aft  body  flexible 
modes  become  more  important.  The  principal  side  bending  moment,,  shear, 
and  torsion  at  Body  Balance  Station  1360  again  results  from  the  Dutch 
roll  response  when  the  yaw  damper  is  off.  All  of  the  frequency  response 
functions  shown  in  Figs.  '27  through  32  are  the  modulus  of  the  complex 
•frequency  responses.  The  load  frequency  response  functions  were  multi¬ 
plied  by  the  gust  power  spectra  which  resulted  in  output  load  power 
spectra.  The  rms  loads  were  obtained  by  taking  the  square  root  of  the 
area  under  the  load  power  spectra.  Typical  variations  of  the  rms  load 
for  the  wing,  fuselage,  and  vertical  tail  are  shown  in  Figs.  33,  35 

and  36,  respectively. 


Determination  of  Critical  Flight  Conditions  and  Critical 
Structure,  Vertical  Analyses 

Initial  efforts  to  arrive  at  design  frequencies  of  exceedance  using  the 
mission  profile  concept  and  to  arrive  at  values  of  ^wld  using  the  design 
envelope  approach  were  based  on  limit  allowable  bending  moments.  The  use 
of  bending  moments  alone  to  estimate  allowable  limit  strength  is  not 
exactly  correct,  because  the  limit  design  stresses  result  from  combined 
loading;  therefore,  properly  combined  shear,  moment,  and  torsion  loads 
should  be  considered.  Since  the  stresses  in  critical  areas  are  more 
sensitive  to  bending  moment,  and  reasonable  assumptions  are  made  as  to 
the  magnitude  and  phasing  of  the  associated  shear  and  torsion,  it  was 
assumed  that  the  critical  structural  areas  and  the  critical  conditions 
could  be  established  by  a  bending  moment  criteria. 

Four  weight  conditions  were  analyzed  in  the  vertical  analyses  for  an 
altitude  of  22,000  ft.  at  the  design  cruise  speed  of  375  kt.  EAS.  The 
resulting  rms  wing  and  tody  bending  moments  for  a  unit  rms  gust  velocity 
were  obtained  and  divided  into  the  incremental  limit  allowable  bending 
moments  to  give  allowable  values  of  These  allowable  <7v/^d  values 

are  plotted  in  Fig.  37  versus  fraction  of  wing  semi-span  and  in  Fig.  38 
versus  Body  Balance  Station.  It  will  be  noted  from  Fig.  37  that  two 
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weight  conditions,  the  fuel  transfer  weight  condition,  186,000  lb.,  and 
the  maximum  flight  weight  condition,  221,600  lb.,  appear  to  be  equally 
critical  inboard  of  33  percent  wing  semi- span.  Since  neither  condition 
was  clearly  critical,  both  weight  conditions  were  further  examined  for 
flight  altitudes  of  5*000,  15,000,  30,000,  and  *<0,000  ft.  The  results 
of  these  altitude  variations  are  given  in  Fig.  39*  wherein  values  of 
A  are  plotted  for  wing  bending  moment  at  Wing  Eta  Stations  .12  and  .33. 
The  corresponding  altitude  variations  of  the  one  fuctor  aeroelastic 
loads  are  given  in  Fig.  40.  The  altitude  conditions  were  analyzed  at 
design  cruise  velocity,  Vg;  thus,  the  true  airspeed  and  Mach  number  in¬ 
creased  with  altitude  until  the  design  cruise  Mach  number  was  reached. 
Above  this  altitude  the  true  airspeed  decreases  for  a  constant  design 
Mach  number  up  to  35*000  ft.  Above  35*000  ft.  the  true  airspeed  remains 
constant  with  increasing  altitude.  The  effect  of  this  variation  of  true 
airspeed  is  apparent  in  the  altitude  trends  for  the  rms  and  one  factor 
flight  bending  moments. 

Values  for  were  calculated  for  each  of  the  two  weight  configura¬ 

tions  mentioned,  previously  for  various  altitudes.  This  information  is 
shown  in  Fig.  4l.  The  discontinuities  in  the  curves  reflect  the  abrupt 
change  in  cruise  velocity  at  24,000  ft.  The  preliminary  boundary  of 
minimum  is  defined  by  Conditions  26,  27,  17,  28  and  21  at  5*000, 

15,000,  22,000,  30,000  and  40,000  ft.,  respectively. 

Speed  variations  to  further  define  the  boundary  of  minimum  aw  V ^  were 
examined  for  altitudes  of  15,000  and  22,000  ft.  The  velocities  chosen 
we re  the  gust  penetration  speed,  VB,  253  kt.  EAS;  an  average  cruise 
speed,  300  kt.  EAS,  and  the  design  dive  speed,  Vp,  445  kt.  EAS  or  Mp, 

.95*  The  fuel  transfer  weight  condition  was  chosen  for  the  22,000  ft. 
speed  variation,  and  the  maximum  flight  weight  was  chosen  for  the 
15*000  ft.  variation.  These  conditions  were  deemed  to  be  the  most  crit¬ 
ical  for  the  respective  altitudes.  The  results  of  the  speed  variations 
are  given  in  Figs.  42  and  43.  The  most  critical  speed  for  the  15*000  ft. 
altitude,  maximum  flight  weight  condition,  appeared  to  be  very  near  the 
design  cruise  speed,  Vq,  375  kt.  EAS  or  472  kt.  TAS;  however,  at  22,000 
ft.  the  critical  speed  was  found  to  be  less  than  Vq.  One  additional 
condition  was  analyzed  for  22,000  ft.  to  pinpoint  the  critical  speed  at 
340  kt.  EAS  or  480  kt.  TAS.  Finally,  center  of  gravity  variations  were 
studied  for  this  condition,  and  the  minimum  (Tvr}(j.  values  were  obtained 
for  a  eg  location  at  .25  MAC,  Condition  24c.  The  resulting  boundary  of 
minimum  Ov/^d  4s  presented  on  Fig.  44.  As  can  be  determined  from  this 
figure,  the'  most  critical  or  least  value  of  <rw  for  Flight  Condition 
24c,  Wing  Eta  Station  .33  is  108.3* 

Minimum  values  of  <rV7>7  &  for  the  fuselage  were  obtained  by  observing  the 
variation  of  ms  fuselage  bending  moments  with  altitude  and  speed.  The 
forward  fuselage  appeared  to  be  more  critical  than  the  aft  fuselage. 

The  one  factor  loads  for  the  forward  body  remained  nearly  the  same  for 
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all  gross  v/eight  conditions;  therefore,  the  variation  in  <rw rj  ^  was 
largely  a  function  of  the  incremental  bending  moments.  The  variation  of 
fuselage  A  values  with  altitude  at  design  cruise  velocity  is  given  in 
Fig.  45(a),.  The  effect  of  speed  variations  on  the  rms  bending  moments 
at  22,000  ft.  is  given  in  Fig.  45(b).  Representing  these  variation  as 
<knd  in  Fig.  46,  it  was  concluded  that  the  vertical  gust  loading  on 
the  fuselage  was  most  critical  at  Body  Balance  Station  480,*  .for  Flight 
Condition  16.  The  minimum  value  of  <TW  rf^  for  this  body  station  was 


112.2. 


The  average  number  of  exceedances  per  hour  of  limit  vertical  bending 
moment  were  calculated  for  the  v/ing  for  Wing  Eta  Stations  .12  and  .33  and 
for  the  fuselage  for  Body  Balance  Stations  480*  and  lloO.  These  results, 
shovm  in  Figs.  47  and  48  respectively,  were  based  on  the  flight  profiles 
describing  the  average  airplane  usage  described  in  Analyses  Conditions. 

The  design  frequency  of  exceedance  for  the  ’.dng  was  determined  as 
3*7(10)-6  exceedances  per  horn.'.  The  corresponding  design  exceedance 
level  for  the  fuselage  is  3*3(lO)"7  exceedances  per  hour.  It  will  be 
observed  that  the  highest  or  most  critical  frequency  of  exceedance  occurs 
at  Wing  Eta  Station  .33,  the  same  station  that  showed  the  minimum  Ow  *7d* 
The  highest  frequency  of  exceedance  for  the  body  also  occurs  at  the  same 
body  balance  station  that  gave  the  minimum  awJ7<i  for  the  fuselage,  Body 
Balance  Station  480.  Therefore,  it  appears  that  these  two  rather  dif¬ 
ferent  approaches  are  consistent  in  arriving  at  the  most  critical  stations 
based  on  the  bending  moment  criteria.  A  comparison  of  the  mission  profile 
and  design  envelope  frequencies  of  exceedance  of  wing  limit  strength  may 
be  of  interest.  By  examining  Fig.  44  herein,  it  is  seen  that  the  crit¬ 
ical  <JwV a  occurs  in  the  altitude  range  22,000  to  27,000  ft.,  and  corres¬ 
ponds  to  an  K(y)/W0  value  of  1.1x10-8.  For  the  critical  condition  the 
bending  moment  zero  crossings  are  1.20  per  sec.,  vrtiich  gives  the  follovj- 
ing  N(y)  value: 

(l. 1x10-8)  (1.20)  (3600)  =  4.8xlO"5  exceedances  per  hour 

As  shown  in  Fig.  47,  the  frequency  of  exceedance  value  as  calculated  by 
the  flight  profile  analysis  is  3.7xlO"6.  The  ratio  of  design  envelope  to 
flight  profile  frequency  of  exceedance  at  the  limit  strength  level,  is 
therefore 


„  3.3.0 

3.7x10-6 

To  extend  the  comparison  to  a  load  basis  let  us  suppose  that  the  flight 
profile  frequency  of  exceedance  were  to  be  applied  on  a  design  envelope 


*This  corresponds  to  Body  Station  540  (See  Fig.  23) 


FIGURE  45.  (a)  VARIATION  OF  FUSELAGE  RMS  VERTICAL  BENDING 

MOMENTS  WITH  ALTITUDE 


A,  FUSELAGE  BENDING  MOMENT  (103  InHb/fp*) 

FIGURE  45.  (b)  VARIATION  OF  FUSELAGE  RMS  VERTICAL  BENDING 

MOMENTS  WITH  SPEED 
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FIGURE  46.  (a)  VARIATION  OF  ALLOWABLE  FUSELAGE  aw  rjd  WITH 

ALTITUDE  BASED  ON  VERTICAL  BENDING  MOMENTS 
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FIGURE  46.  (b)  VARIATION  OF  ALLOWABLE  FUSELAGE  <rwj;d  WITH 
SPEED  BASED  ON  VERTICAL  BENDING  MOMENTS 
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basis  -  for  example,  in  order  to  provide  for  operation  of  the  airplane 
100$  of  the  time  at  its  most  critical  flight  condition.  The  N(y)/N0 
value  for  the  critical  design  envelope  condition  using  the  flight  pro¬ 
file  N(y)  with  the  design  envelope  N0  is: 

^  “  (3-7x10-6)  (_!_)  (_!_)  „  8.6x10-1° 

At  an  altitude  of  22,000  ft.,  the  corresponding  Ov1}  from  Fig.  47  is 
138  ft.  per  sec.  But  -strength  has  been  provided  only  for  <rwrj  ^  =  108.3 
ft.  per  sec.  As  a  result,  the  incremental  loads  would  have  to  be  in¬ 
creased  in  the  ratio  I38/IO8.3  =  1.27.  The  corresponding  ratio  of  ^ 
increase  of  net  loads,  assuming  a  one-factor  bending  moment  of  8.00x10° 
in. lb.  and  a  limit  allowable  bending  moment  of  22.25x10°  in. lb.,  is 

(22.25  -  8.00)  (1.27)  +  8.00  _  . 

22.25  “ 


Determination  of  Critical  Flight  Conditions  and  Critical  Structure, 
Lateral  Analyses 

Previous  lateral  analyses  of  707/720  type  airplanes  have  indicated  that 
the  magnitudes  of  the  gust  loads  for  flight  in  continuous  turbulence  is 
largely  dependent  on  the  Dutch  roll  stability  of  the  airplane.  Naturally 
the  loads  are  influenced  by  the  damping  in  the  Dutch  roll  mode,  and  can 
reach  large  amplitudes  as  the  damping  decreases.  There  are  regions 
within  the  flight  regime  where  the  Dutch-  roll  damping  decreases  to  the 
extent  that  the  Dutch  roll  response  would  reach  significant  amplitudes 
if  the  controls  were  locked.  However,  the  pilot  would  apply  corrective 
control  inputs  either  manually  or  automatically  through  the  yaw  damper 
system  before  such  a  situation  would  develop.  Since  this  corrective 
action  is  at  the  discretion  of  the  pilot,  it  is  beyond  the  scope  of  this 
study  to  establish  a  threshold  for  corrective  control. 

Therefore,  the  lateral  analyses  for  the  design  envelope  approach  were 
conducted  for  flight  conditions  defined  by  the  design  cruise  velocity, 

Vc,  at  a  variety  of  altitudes.  An  airplane  configuration  was  chosen  to 
give  the  greatest  Dutch  roll  response.  This  configuration  exhibits  the 
maximum  airplane  yaw  inertia,  a  high  gross  weight  and  a  large  yaw-roll 
product  of  inertia,  all  of  which  are  associated  with  reduced  Dutch  roll 
stability.  The  rms  fin,  and  aft  fuselage  bending  moments  were  divided 
into  the  limit  allowable  bending  moments  for  the  fin  and  aft  fuselage, 
respectively,  to  give  values  of  OyVd’  ^hese  <?v^d  values  were  plotted 
against  Fin  Elastic  Axis  Station  and  Body  Balance  Station,  and  are  shown 
in  Figs.  49  and  50.  The  minimum  value  of  OV^d,  6L*9»  for  the  fin  occurs 
at  Fin  Elastic  Axis  Station  158.  The  most  critical  aft  body  balance 
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station  is  1080  where  the  <rv  *?  value  was  9Jt*  3 •  The  critical  aft  fuse¬ 
lage  station  from  the  lateral  analysis  standpoint  is  very  nearly  the  same 
as  critical  aft  fuselage  station  determined  from  the  vertical  analyses-; 
Resultant  values  of  <rv  17  ^  obtained  for  other  altitude  points  on  the  de¬ 
sign  cruise  speed  envelope  are  shown  in  Figs.  51  and  52.  The  critical 
altitude  for  both  the  aft  fuselage  and  fin  is  23,500  ft.  All  of  the 
above  analyses  conditions  were  repeated  with  the  electronic  yaw  damper 
engaged.  It  will  be  noted,  from  the  results  shown  in  Figs.  51  and  52  that 
the  yaw  damper  is  very  effective  in  reducing  the  peak  amplitude  of  the 
Dutch  roll  response,  and  therefore  provides  a  large  reduction,  approxi¬ 
mately  40  to  60  percent,  in  the  fin  and  fuselage  loads.  The  degree  to 
which  the  ya w  damper  is.  effective  in  reducing  the  loads  is  largely 
dependent  upon  the  effective  Dutch  roll  damping  it  can  provide. 

The  number  of  exceedances  of  limit  design  bending  moment  per  average 
flight  hour  were  calculated  for  the  vertical  tail  Elastic  Axis  Stations 
6?  and  158  and  Aft  Fuselage  Balance  Stations  1020  and  1335*  These  values, 
obtained  both  with  the  yaw  damper  engaged  and  disengaged,  are  shown  in 
Figs.  53>  54,  55,  and  56.  The  most  critical  design  value  for  the  verti¬ 
cal  tail,  1.05(10)“^  exceedances  per  hour,  occurs  at  Elastic  Axis  Station 
158  with  the  yaw  damper  disengaged.  When  the  yaw  damper  is  operating 
full  time  the  value  is  increased  to  4.0(l0)~6  exceedances  per  hour,  thus 
showing  the  capability  of  yaw  damper  system  to  reduce  the  level  of  the 
lateral  dynamic  loads.  The  design  frequency  of  exceedance  for  the  aft 
fuselage  occurs  at  Body  Balance  Station  1020  and  with  the  yaw  damper 
disengaged  is  2,0('10)“'  exceedances  per.  hour.  The  corresponding  values 
with  the  yaw  damper  engaged  is  1.50(lO)"9  exceedances  per  hour. 

The  critical  or  design  values  of  <Tvrj  ^  based  on  bending  moment  from  the 
lateral  analyses  occur  at  the  same  vertical  tail  .station  and  aft  body 
station  as  the  critical  design  frequency  of  exceedance  obtained  by  the 
flight  profile  approach. 

A  summary  of  the  exceedances  of  limit  design  load,  as  obtained  from  the 
Flight  Profile  Analyses  using  the  bending  moment  criteria,  are  given  in 
Fig.  57*  Also  shown  in  this  figure,  expressed  as  exceedance  ratios,  are 

17  d )  N(^v^d^ 

the  design  values  of  - r= -  and  - ^ -  as  determined  from  the  Design 

*’0  ^0 

Envelope  Analyses. 


Results  -  Incremental  CG  Acceleration 

The  rms  values,  A,  of  incremental  center  of  gravity  acceleration  due  to 
unit  ms  gust  velocity  are  shown  in  Fig.  58  for  gross  weights  of  221,600 
lb.  and  186,000  lb.  These  values  were  obtained  from  the  dynamic  loads 
solution  for  points  on  the  design  cruise  speed  boundary  of  375  kt.  EAS 
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FIGURE  53.  EXCEEDANCES  OF  FIN  BENDING  MOMENT,  YAW 
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FIGURE  54.  EXCEEDANCES  OF  FIN  BENDING  MOMENT,  YAW 
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FIGURE  56.  EXCEEDANCES  OF  FUSELAGE  LATERAL  BENDING 
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FIGURE  57.  SUMMARY  OF  DESIGN  VALUES  OF  EXCEEDANCE  RATIOS 
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to  Mach  .90.  As  a  matter  of  interest,  the  A  values  for  186,000  lb.  gross 
weight  were  also  calculated  using  the  one  dimensional  expression  developed 
by  Y.  C.  Fung(To)  and  given  below: 

pvToLa 
'^An  =  2(w/s) 

where, 

p  atmospheric  density  at  altitude, 

slugs  per  cu.  ft. 

true  airspeed,  'ft.  per  sec. 

flexible  airplane  life  curve  slope, 
a  per  rad. 

W  airplane  gross  weight,  lb. 

S  airplane  wing  area,  ft. 

gust  response  factor  based  upon  the  isotropic 
turbulence  spectrum  and  a  scale  of  turbulence 
factor  of  2500  ft. 

The  A^  values  determined  by  the  one-dimensional  equation  are  about  ten 
percent  higher  than  those  obtained  from  the  full  dynamic  solution.  This 
difference  can  be  attributed  to  flexibility  and  aerodynamic  effects  which 
are  more  precisely  accounted  for  in  the  dynamic  solution. 

Values  for  the  frequency  of  exceedance  per  hour  for  incremental  center 
of  gravity  acceleration  are  shown  in  Fig.  5Sh  These  values  were  deter¬ 
mined  from  flight  profile  analysis.  The  eg  acceleration  exceedances  for 
the  707-300  series  airplanes  are  also  shown  in  Fig.  59*  These  data  were 
collected  by  NASA  from  actual  airline  operations^),  At  the  time  of  this 
report,  sufficient  acceleration  data  was  not  available  for  the  720B  air¬ 
plane,  so  an  attempt  was  made  to  adjust  the  707-300  data  to  represent 
720B  usage.  This  adjusted  data  can  be  compared  to  the  720B  analytical 
data  in  Fig.  59*  This  comparison  indicates  that  the  predicted  eg  accel¬ 
eration  exceedances  are  somewhat  higher  than  the  adjusted  usage  values. 

The  disagreement  increases  as  the  acceleration  level  increases.  This 
would  indicate  that  the  turbulence  intensity  parameters,  Pp,  P2,  b]_,  and 
b2  should  be  readjusted. 


Ill 
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Boeing  Analysis  Results  for  Lockheed  Model  7^9 

A  comparative  analysis  was  conducted  by  The  Boeing  Company  on  the  Lock¬ 
heed  Model  7^9  airplane.  The  Lockheed-California  Company  furnished  the 
basic  airplane  data  in  the  form  of  lumped  mass  distribution,  bending  and 
torsional  stiffness  distributions,  wing  lift  distribution,  and  basic 
airplane  lift  and  moment  data.  From  this  description  the  analysis  was 
condir  'd  by  Boeing  using  the  analysis  procedures  used  for  the  720B 
study.  The  resuitr*  pf  this  analysis  are  expressed' as  A,.,  the  rms  load 
due  -to  a  one  It.  per  sec.  rms  gust  velocity,,  and  N0,  the  number  of  zero 
crossings  per  second  for  wing  shear,  bending  moment,  and  torsion.  A 
comparison  of  these  values  with  those  obtained  by  Lockheed  are  shown  in 
Figs.  60  and  6l.  Values.  for  A  and  N0  for  eg  incremental  acceleration 
are  also  listed.  The  agreement  between  the  Boeing  and  Lockheed  A  values 
is  generally  not  as  good  as  might  be.  desired.  Possible  reasons  for  this 
disagreement  are  explored  in  Reference  18.  The  agreement  of  the  N0 
values  is  seen  to  be- much  more  satisfactory. 


Parametric  Variations  -  Vertical  'Analyses 

The  results  of  the  parametric  variations  conducted  for  the  vertical  anal¬ 
yses  are  shown  in  Fig.  §2  as  they  affected  the  rms  wing  bending  moment 
at  the  critical  wing  station.  Fig.  63  shows  how  these  parametric  var¬ 
iations  affect  the  airplane  incremental  eg  acceleration.  The  trends, 
'exhibited  by  bending  moment  and  eg  acceleration  are  similar  for  all 
variations  except  for  the  variation  of  wing  torsional  stiffness.  In 
this  case,  a  reduction  in  wing  torsional  stiffness  of  20  percent  re¬ 
sulted  in  a  slight,  less  than  one  percent,  increase  in  eg  acceleration 
while  the  rms  wing  bending  moment  actually  decreased  by  one  percent. 

The  wing  rms  bending  moment  and  incremental  eg  acceleration  are  quite 
insensitive  to  variations  in  wing  torsional  stiffness. 

When  gradual  penetration  effects  were  neglected,  the  resulting  rms  wing 
bending  moment  and  the  airplane  incremental.  eg-  acceleration  were  reduced 
by  10-12  percent.  This  implies  that  results  obtained  from  analyses,  in 
which  gradual  penetration  was  neglected  may  be  unconservative.  The  rms 
wing  bending  moments  and  incremental  eg  acceleration  were  also  affected 
by  the  inclusion  of  structural  damping.  Including  structural  damping  in 
the  amount  of  g  =  .03  reduced  the  rms  wing  bending  moment  and  eg  accel¬ 

eration  levels  from  5  to  8  percent.  Moving  the  airplane-  center  of 
gravity  aft  by  9  percent  of  the  MAC  increased  the  wing  bending  moment 
and  eg  acceleration  2  to  3  percent. 

The  rms  wing  bending  moment  level  was  reduced  about  8  percent  when  the 
wing  bending  stiffness  was  decreased  to  80  percent  of  its  nominal  value. 


INCLUDING 


The  general  conclusions  that  cun  be  drawn  from  the  vertical  analyses 
par  rone trie  variations  is  that  the  dynamic  analyses  should  include  gradual 
penetration  and  structural  damping  effects. 


Parametric  Variations  -  Lateral  Analyses 


The  results  of  the  parametric  variations  conducted  for  the  lateral  anal¬ 
ysis  are  given  in  terms  of  vertical  tail  bending  moment  at  .Fin  Elastic 
Axis  Station  158.  These  results  are  shown  in  Figs.  64  through  66(c). 
Variations  in  A  for  fin  shear,  bending  moment,  and  torsion  are  shown  in 
Fig.  64  for  changes  in  airplane  gross  weight.  The  gross  weight  changes 
were  made  by  adding  wing  fuel.  The  addition  of  wing  fuel  causes  a 
corresponding  increase  in  airplane  yaw  and  roll  inertia.  The  trends  in 
the  rms  loads,  therefore,  reflect  the  influence  of  these  inertia  changes 
as  well  as  the  weight  changes.  These  results  indicate  a  11  percent  in¬ 
crease  in  the  rms  fin  loads  for  a  25  percent  increase  in  gross  weight. 


The  variation  in  A  for  fin  loads,  resulting  from  changes  in  altitude, 
with  gross  weight,  Mach  number  and  true  airspeed  constant  at  475  kti, 
are  shown  in  Fig..  65.  The  effects  of  increasing  altitude  under  these 
circumstances  is  to  decrease  the  dynamic  pressure  and  decrease  Dutch  roll 
damping.  The  reduction  of  Dutch  roll  damping  would  tend  to*  increase  the 
resulting  loads;  however,  for  this  airspeed,  the  reduction  of  dynamic 
pressure  more  than  offsets  the  effect  of  the  damping  loss,  and  the  end 
result  was  an  overall  decrease  in  the  loads.  Previous  studies  on  similar 
airplanes  indicate  that  increasing  the  altitude  at  reduced  airspeeds  can 
increase  the  loads  significantly  due  to  a  loss  of  Dutch  roll  stability. 


The  curves  in  Figs.  66(a)  through  66(c)  show  the  effects  on  rms  fin  loads 
of  variations  in  the  yaw  damping  derivative,  Cn.  •  the  yaw-roll 


coupling  derivative  y  Cn. 


'b^/2V' 


■'nb^/2V> 

and  the  airplane  product  of  inertia, 


xy 


As  was  mentioned  previously  in  Analysis,  the  lateral  gust  loads  are 
closely  associated  with  the  Dutch  roll  damping.  When  the  Dutch  roll 
damping  is  large,  variation  in  the  damping  derivatives  and  airplane  pro¬ 
duct  of  inertia  have  only  a  small  effect  on  the  lateral  dynamic  loads; 
however,  when  the  amount  of  inherent,  '"'itch  roll  damping  is  low,  these 
same  parameters  can  have  a  tremendous  effect  on  the  dynamic  stability  and 
the  resultant  dynamic  loads. 


Statistical  Correlation  Between  Wing  Loads 


Prior  to  conducting  the  detailed  combined  stress  analysis  it  was  decided 
to  obtain  the  correlations  between  the  wing  shear,  bending  moment  and 
torsion  loads.  These  correlation  coefficients  were  obtained  for  each 
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A,  FIN  TORSION  (10  0  in. -lb  fpj) 


FIGURE  66  (c).  FIN  TORSION  RESULTS  OF  PARAMETRIC  VARIATIONS, 

LATERAL  ANALYSES 


set  of  load  quantities  vising  cross  power  spectral  density  equations 
•derived  in  Appendix  B.  The  results  are  , plotted  versus  Wing  Eta  Station 
in  Fig.  67. 

Outboard  of  the  outboard  nacelle  the  wing  shear  and  bending  moment  ai?e 
almost  exactly  correlated  while  the  torsion  and  bending  moment  axe 
correlated  very  poorly.  Hiis  indicates  that  shear  and  bending  are  very 
nearly  in  phase  but  that  torsion  is  very  nearly  independent  of  either 
bending  moment  or  shear;  Inboard  of  the  outboard  nacelle,  the  effects 
of  the  nacelle  are  readily  apparent.  The  correlation  between  shear  and 
bending  moment  is  poorer  due  to  shear  loads  arising  from  the  outboard 
nacelle.  However,  the  correlation  between  shear  and  torsion  is  improved, 
due  to  the  nacelle  inertia  loads.  A  somewhat  similar  tendency  is  shown 
at  the  inboard  nacelle.  The  correlation  between  the  three  load  quan¬ 
tities  improves  inboard  of  the  inboard  nacelle.  This  is  a  result  of 
transferring  shear  and  bending  moment  into  torsion  as  the  elastic  axis 
sweep  angle  changes  near  the  wing  root. 


Results  -  RMS  Stress  Analysis.. 

The  individual  structural  elements  for  the  detailed  stress  analysis  were 
selected  at  locations  around  the  critical  wing,  body  and  fin  stations  to 
give  the  best  strength  evaluation  possible.  Included  among  these  ele¬ 
ments  are  those  which  demonstrated  the  least  margin  of  safety  as  deter¬ 
mined  by  the  static  stress  analysis.  (11>12,13)  The  numbering  system 
used  for  all  the  elements  is  the  same  as  that  used  in  the  stress  anal¬ 
ysis  referenced  above.  In  Fig.  68  is  a  sketch  identifying  the  structural 
elements  of  the  wing  box  at  Wing  Eta  Station  .33.  The  elements,  for  which 
design  values  of  and  number  of  exceedances  were  calculated  are 

indicated  by  the  solid  circles.  The  structural  element  array  and  the 
elements  analyzed  for  Body  Stations*  540,  1040,  and  1360  are  show,  in 
Figs.  69(a),  69(b),  and  69(c),  respectively.  The  same  elements  were 
analyzed  in  both  the  lateral  and  vertical  analyses  with  the  exception 
of  Forward  Body  Station  540  which  was  analyzed  for  the.  vertical  analyses 
only.  The  detailed  stresses  for  the  critical  Vertical  Tail  Elastic  Axis 
Station  158  were  calculated  for  the  elements  shown  in  Fig.  70.  The 
stress  analysis  for  the  fin  was  conducted  assuming  the  elements  on  the 
left  and  right  side  were  symmetrical-  about- the  fin  chord  line. 

The  rms  loads  were  transformed  into  axial  and  shear  stresses  using  the 
stress  coefficients  obtained  from  the  conventional  stress  analysis.  The 
rms  stress  quantities  for  the  structural  elements  at  Wing  Eta  Station 
•33;  Condition  24c,  are  shown  in  Tables  10(a)  and-10(b).  The  rms  §hear, 


*Corresponds  to  Body  Balance  Stations  480,  1000,  and  1320  (see 'Fig.  21) 


FIGURE  68.  ARRAY  OF  STRUCTURAL  ELEMENTS,  WING  ETA 

STATION  .  33 


M 


NO  SCALE  pr\  STRUCTURAL  ELEMENTS 
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FIGURE  69  (a).  ARRAY  OF  STRUCTURAL  ELEMENTS,  BODY 

STATION  540 
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FIGURE  69  (b).  ARRAY  OF  STRUCTURAL  ELEMENTS,  BODY 

STATION  1040 
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ARRAY  OF  STRUCTURAL  ELEMENTS,  BODY 
STATION  1360 


AIRPLANE 
PLAN  VIEW 


STRUCTURAL 

ELEMENTS 

ANALYZED 


NO  SCALE 

FIGURE  70.  ARRAY  OF  STRUCTURAL  ELEMENTS,  FIN  ELASTIC  AXIS 

STATION  158 


'iiliisiiiiiiiiffiiiiiii-ifffiisiiiiisiiiifis 


•  •  -  ' 

s  *5  lif f siif Issiissif isssisisilf liiiiisiiiisiii 


#•  ••#••••••••••••♦••  •-•  •  ••  •.#  ••••••••••  ••••••  •  • 


■m  *«•#*»«•#»# •  *«•»«**  f  •  «*«*»*  •«««*#*  •  ••  •£#££># 


*  «n»>«  »+»###•>•> 

I»  ••••  •••  ••%•  •  •••  •  •  •  •  ••••  •  »  *  •#•••••••  ••••••••• 


10b.  CORRELATION  COEFFICIENTS  AND  STREJ 
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skin,  and  segment  stress  levels,. A,  are  given  in  Table  10 (a)  with  their 
associated  number  of  zero  crossings.  Table  10(b)  contains  the  correla¬ 
tion  coefficients  between  the  shear  and  axial  stresses  and  the  respective 
stress  rates.  Similar  values  were  calculated  for  all  the  structural 
elements  indicated  in  Figs.  68  through  70. 


Results  -  Statistical  Analysis 


Following  the  procedure  outlined  in  the  Analysis  section,  the  rms  stresses, 
correlation  coefficients,  and  stress  rates  were  used  to  determine  the 
probability  of  exceeding  limit  strength,  P(MS  <  0,  <XW)  ,  and  the  number  of 
times  per  hour  that  -the  limit  strength  envelope  is  penetrated.  The  re¬ 
sultant  probabilities,  P(MS  <0,  tfw) ,  are  given  in  Table  11  for  Structural 
Element  9>  Condition  24c.  The  probability,  P(MS<0,  <jrw),  for  all  the 
structural  elements  analyzed  by  the  design  envelope  approach  were  plotted 
and  are  shown  in  Figs.  7l(a)  through  7^(d).  The  most  critical  element 
for  each  profile  condition,  i.e.,  the  element  which  demonstrates  the 
highest  probability  of  exceeding  limit  strength,  was  selected  using  the 
techniques  described  on  Page  58.  The  value^  resulting  from  the  multipli¬ 
cation  of  these  critical  P(MS<0,  <rw)  with  2  (<TW)  are  shown  in  Figs.  75 
through  78.  The  maximum  values  ofAthese  overall  probability  curves  were 
used  with  the  curves  of  P(An  >An^)f (Oy)  in  Figs.  79  and  80  to  determine 
the  design  incremental  eg  acceleration.  The  design  levels  of  ^7  a  and 
Ov  7a  were  then  calculated  and  appear  in  Table  12.  Corresponding  exceed¬ 
ance  ratios  as  determined  by  Eq.  (15 )  are  also  given  in  Table  12. 

The  altitude  variation  of  <rw  determined  by  the  joint  probability 
method  is  shown  in  Fig.  8l  for  Wing  Eta  Station  .33  and  for  the  vertical 
-tail  E.  A.  Station  158.  Similar  information  is  given  in  Fig.  82. 

The  altitude  trends  obtained  by  the  joint  probability  approach  indicate 
that  the  critical  altitude  of  23,p00  ft.  for  the  lateral  analyses  and 
22,000  ft.  for  the  vertical  analyses  as-  chosen  by  the  bending  moment 
criteria  were  correct. 

The  variation  of  *w7d  around  the  wing  structural  box  at  Wing  Eta 
Station  .33,  Condition  24c,  as  determined  by  the  joint  probability 
method,  is  -shown  in  Fig.  83.  The  static  margin  of  safety  for  each  ele¬ 
ment  as  obtained  by  the  static  stress  analysis  (12 )  for  a  high  bending 
moment  condition  is  also  shown  in  Fig.  83«  The  comparison  indicates 
that  the  strength  capability  as  determined  by  the  statistical  approach 
e;;hibits  the  same  trend  as  the  static  margin  of  safety.  There  is,  how¬ 
ever,  a  difference  in  which  element  is  most  critical.  Element  9  and  the 
lower  surface  has  the  lowest  margin  of  safety  while  Elements  117  and  122 
on  the  upper  surface  display  the  lowest  <rw7  d*  The  significance  of  this 
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TABLE  11.  VALUES  OF  P(MS'O)  f  (cw) 
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PROBABILITY  OF  EXCEEDING  ZERO  MARGIN.  P(MS<0. 


RMS  GUST  VELOCITY  (ft/wc) 

FIGURE  71  (b).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN 

CONDITION  24 


PROBABILITY  OF  EXCEEDING  ZERO  MARGIN,  P(MS<0,  <xj 


V  RMS  GUST  VELOCITY  (ft/»*c) 

FIGURE  71  (c).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  24C 
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PROBABILITY  OF  EXCEEDING  ZERO  MARGIN,  P(MS<0,  a  ) 


PROBABILITY  OF  EXCEEDING  ZERO  MARGIN,  P(MS<0 ,o 
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FIGURE  72  (a),  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  16 


PROBABILITY  OF  EXCEEDING  ZERO  MARGIN,  P(MS  <  0, 


FIGURE  72  (b).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  16 
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PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  P(MS<0, 


FIGURE  73  (a).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN 

CONDITION  3B 
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FIGURE  73  (b).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  4B 


PROBABILITY  OF  EXCEEDING  ZERO  MARGIN.  P(MS<0 


FIGURE  73  (c).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN 

CONDITION  5B 


PROBABILITY  OF  EXCEEDING  ZERO  MARGIN,  P(MS<0.  <Ty) 
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FIGURE  73  (d).'  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  3BYD 
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FIN  EA  STATION  158 
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FIGURE  73  (e).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN 

CONDITION  4BYD 


PROBABILITY  OF  EXCEEDING  ZERO  MARGIN,  P^S^O^ 


FIGURE  73  (f).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  5BYD 
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FIGURE  74  (a).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  -  » 

CONDITION  4B  i 
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FIGURE  74  (b).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  4B 
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FIGURE  74  (c).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  4BYD 


PROBABILITY  OF  EXCEEDING  ZERO  MARGIN,  P(MS  <0 


FIGURE  74  (d).  PROBABILITY  OF  EXCEEDING  ZERO  MARGIN  - 

CONDITION  4BYD 


WING  rj  STATION  .33 


(rw,  RMS  GUST  VELOCITY  (fps) 

FIGURE  75.  OVERALL  PROBABILITY  OF  FAILURE  FOR  CRITICAL 
ELEMENTS-  WING;  VERTICAL  ANALYSIS 
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Values  of  Px  f  (dL)  plotted  In  this  figure  ore  based 
on  one  dimenslonol  value  of  P  ( MS  £  5)  vs  <r*.  The 
reason  Is  that  the  correlation  coefficients  are  very 
near  unity  and  the  predominance  of  one  stress  over  the 
other  has  led  to  inconsistent  numerical  integration 
of  the Joint  probability  function*  It  is  believed  that 
the  ono  dimensional  value  is  sufficiently  accurate  to 
represent  aliplane  strength. 


body  Vertical  analysis 

CONDITION  16 
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<rWf  RMS  GUST  VELOCITY  (fpsj 

^  ;  FIGURE  76.  OVERALL  PROBABILITY  OF  FAILURE  FOR  CRITICAL 

ELEMENTS- BODY,  VERTICAL  ANALYSIS 
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P(MSs  0)  f( 


BODY  LATERAL  ANALYSIS 


Valuet  of  Px  t  (<rv)  plotted  In  this  figure  are  based 
on  one  dimensional  value  of  P  (MS  £  0)  vs  <rv.  The 
reason  is  that  the  correlation  coefficients  are  very 
pear  unity  and  the  predominance  of  one  stress  over  the 
other  has  led  to  inconsistent  numerical  integration 
of  the  joint  probability  function.  It  is  believed  that 
the  one  dimensional  value- is  sufficiently  occurate  to 
represent  otrptane  strength. 


VtLUUTY  (fps) 

FIGURE  77.  OVERALL  PROBABILITY  OF  FAILURE  FOR  CRITICAL 
ELEMENTS- BODY,  LATERAL  ANALYSIS  UKITI°AL 
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P(MS  *0)f(a  ) 


FIN  EA  STATION  158 


<TV,  RMS  GUST  VELOCITY  (fp*) 


FIGURE  78.  OVERALL  PROBABILITY  OF  FAILURE  FOR  CRITICAL 
ELEMENTS- VERTICAL  TAIL,  LATERAL  ANALYSIS 
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FIGURE  81. 


VARIATION  OF  COMBINED  STRESS  PROBABILITY  *wi»d 
WITH  ALTITUDE 
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ALTITUDE,  1000  FT 


An. 

“A  =  "vld 

FIGURE  82.  VARIATION  OF  COMBINED  STRESS  PROBABILITY  *v  >>d 

WITH  ALTITUDE 
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FIGURE  83.  COMPARISON  OF  STATIC  MARGIN  OF  SAFETY  WITH 
v*  -  -  WING  ETA  STATION  .  33 


-difference  is  that  an  allowable  bending  moment,  for  the  section  would 
probably  be  determined’ by  Element  9  while  actually  under  a  random  type 
of  loading  Elements  117  and  122  may  be  the  most  critical. 

The  evaluation  of  strength  capability  following,  the  Flight  Profile 
Approach  results  in  number  of  exceedances  per  hour  of  the  limit  strength 
envelope.  Representative  values  of  the  number  of  exceedances  per  hour, 

G,  for  structural  Element  110  are  given  in  Tables  13  through  25  for  all 
the  flight  profile  conditions.  The  G  values  for  the  remaining  elements 
and.  conditions  appears  in  Volume  II,  Da.ta  Report.  The  total  number  of 
exceedances  of  the  limit  design  envelope  per"  average  flight  hour,  Np, 
were  obtained -by  multiplying ■ the  G  values  by  the  appropriate  time 
factors.  These  values  are  summarized  in  Tables  26(a),  26(b)  and  2 6(c). 

A  summary  of  the  exceedances  per  hour  obtained  for  the  Flight  Profile 
Analysis.-  using  the  joint  probability  approach  are  given  in  Fig.  84.  Also 
shown  in  this  figure  are  the  design  values  of  tfw^d  and  Oy^d  expressed 

N(<rwrfd)  N(<rvi?d) 

as  exceedance  ratios  — ^ —  and  — —  * 
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TABLE  13.  EXCEEDANCES  PER  HOUR -  CONDITION  2 

7208  CONDITION  2  VERTICAL  COST  ANALYSIS 
MING  ETA  STATION  •  0.31  STRUCTURAL  EL2N2NT  110 
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TABLE  14.  EXCEEDANCES  PER  HOUR  -  CONDITION  4 
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TABLE  15*  EXCEEDANCES  PER  HOUR  -  CONDITION  5 

not  COMO  IT  ION  S  VERTICAL  COST  ANALYSIS 
N1NCETA  STATION  *0.39  STRUCTURAL  2LFMENT  HO 
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0.29073182-07 

99. 

0.18776972-01 

6.27892292-10 

0.926970.2-08 

80. 

0.27038722-01 

0.21989262-11 

0.49662472-09 

89. 

0.99219082-01 

0.11177071-12 

0.99007992-10 

70. 

0.92707882-03 

0.66392012-14 

0.9911899E-U 

79. 

0.6898677E-01 

0.26146192-19 

0.19891792-12 

80. 

0.81890002-01 

0.83741102-17 

0.7691299Erl4 
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TABLE  17.  EXCEEDANCES  PER  HOUR  -  CONDITION  7 


MOB  C0N0IT10N  7  VERTICAL  GUST 

ANALYSIS 

MING  ETA 

STATION  «  0.33  STRUCTURAL  ELEMENT  110 

TURBULENCE 

PARAMETERS 

PI  •  0.016S00 

P2  •  0.000090 

61  •  3,050000 

82  *11,400000 

1G  STRESSES  IPSO 

SHEAR  SKIN 

SEGMENT 

- 

-2274.0  -11740.0 

-11220.0 

V  TRUE  IPT/SEC 1  *010.7 

G  BAR  »  0. 

24IRSS7E-05  PER  HOUR 

-RNS 

NC 

f  ISIGM) 

NC  fISIGMT 

GUST  VELOCITY 

X  1000V 

If T/SEC 1 

IPER  FT  1 

10. 

0. 

0.2427B38E-04 

0. 

IS. 

0. Ill  9441 F-20 

0.2674668E-OS 

0.4369? 1 7F-20 

20. 

0.1 36TS00F-1? 

0.13SIS4SF.-0S 

0.26970 34F-12 

22. 

0.8  729709E-1 1 

0.9TRS?  T4F-06 

0. 1 246S3RE- 10 

24. 

0.2060932E-09 

0.686834TF-06 

0.206S612E-09 

26. 

0.24141 29E-08 

0.4674R16F-06 

0.1646R64E-08 

26. 

0. 1 701502E-0T 

O.3OH4300F-O6 

O.T66O014E-O0 

SO. 

0. R223965E-0T 

0.1974640F-06 

0.2369746E-07 

3?. 

0.29R6331F-06 

0.  1  224440F -06 

0.S3403S1 E-07 

34. 

0.8696310E-06 

O.T374602F-OT 

0.935B624E-07 

36. 

0.2129936E-05 

O.43O3407F-O7 

0. 1 337580F-06 

3A. 

0.454S960E-05 

0.243SI 77E-07 

0.161S433E-06 

40. 

0. 048481 6F -OS 

O.13362O7F-07 

0.1693429F-06 

4$. 

O.3OS4043F-O4 

O.26O4640F-O8 

0. t 161 1 07F-06 

SO. 

0. 7SI 2609E-04 

0.4 1  PR  745F-09 

0.4S92059E-07 

SS. 

0.1462309F-03 

0.5S47420E-IO 

0.  1 105092F.-O7 

60. 

0. 242  T509f -03 

0.600301 2F-1 1 

O.2IS4026E-O8 

6S. 

0.  36O20O2F-O3 

0.5493U34E-1? 

0.200  79  74F-09 

TO. 

0.4931 336F-03 

0.409?  390F - 1  3 

0.2944931  E-  TO1 

TS. 

0.43SR29«E-03 

0.251S31 IF-14 

0.2333609F-II 

no. 

0.  7030096E-O3 

0.12  TS44 SF- t  S 

0. I4S090OF- I 7 
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TABLE  ,18.  EXCEEDANCES  PER  HOUR  -  CONDITION  7  (31,000  FT. ) 


>t?oo  conoition  t  vfotical  oust 

analysis 

WING  ETA  STATION  •  0.31  STOIfC  TUOAL  Elf  HE  NT  1 10 

ALTITUDE  31,000  FT. 

TU«0UtENCF  OAOANETFOS 

01  •  O.OISOOO 

0?  •  0.000000 

01  •  3.000000 

0?  «11.7SO000 

10  STOESf.FS 

IF  SI  1 

SHfAO  SKIN 

SFGMFNT 

-7.7  IA.0  -IITAO.O 

-11770.0 

V  T0UE  IFT/SFC1  *010. T 

0  OAF  •  0.20SMA9r-0S  PF*  MOWO 

•NS 

NC 

F  ISir.Ml 

NC  FISIGM) 

OUST  VELOCITY 

1  IOOOV* 

iFT/srct 

IPfO  Ml 

10. 

0. 

0.I9T776TF-0A 

0. 

is. 

0.111  'IAA  l  E  ”  70 

0.  2A39078E -OS 

0.A31099TE-70 

70. 

0.  1  3  A  TS  BOf  - 1  ? 

0.  1 31 AA  31  E-OS 

0.76731ASE-12 

77. 

0.BT?0T0OF-U 

0.9A3776TE-0A 

0.  1 701 SAOF-  1 0  .. 

7A. 

0. 70609 1?F -09 

0.6SST961F-0A 

0.I9T77ASE-09 

7A. 

0.7A1 AI79F-00 

0.AA1 TA9SF-06 

O.ISSA70AE-00 

?0. 

0.TT01S02F-0T 

0.70033AAE-06 

0. T1S91STF-00 

10. 

0.0??  39ASF--0T 

O.10733S0F-OA 

0.2I00I9AE-0T 

Vi 

0.  79BA33I E -06 

0.11171 TSF-00 

O.A060AA3F-OT 

3*. 

O.OA96330F-OA 

0. A637002E -OT 

O.OAI6IATE-OT 

30. 

0.  ?  1  799  IAF-OS 

0. 301 ASAIE-OT 

0.1 I0SA09F-0A 

3*. 

0.SSAS9A0C-0S 

0.717STSOF-0T 

0.1A10I TIE-06 

AO. 

O.FAflAOJAC-OS 

0.11ATT03F-0T 

0. IASAA37E-OA 

AS. 

0. 30SA8A3E-0A 

0.?l A170AE-00 

0.9SASA3SF-0T 

so. 

O.TS17609F-0A 

0.  3?  TOTOIF-09 

0. 3S9A30IF-OT 

ss. 

O.|AA?109r-O3 

0.A170A3AF-I0 

O.OT9?SATf-O0 

PO. 

0.2A7  TSmF-01 

0.A7S0003F-U 

O.1SOSS3SF-O0 

AS. 

o.tao?fo?e-oi 

0.359B03SF-I? 

0.1 09| AAOF -09 

TO. 

O.A931 33AF-03 

0.7S00039F- 1 3 

0.1T99OS1F-1O 

TS. 

O.A  3S8798E -03 

0. 1A7STAAF -IA 

O.13720OAE-I1 

00. 

0. T83K09AF-03 

0.AAT3AAI F-IP 

0.T613TA0F-13 
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TABLE  19.  EXCEEDANCES  PER  HOUR  -  CONDITION  8 


770*  coNomnN  A  vK»»if.*i  oiim 

ANA  I  VS  RS 

WINf.  fi*  MAIION  .  0.11 

MPUCNJAAi  f  1  f  Nf  N1  110 

VUMIUrNCF  PAAAMfTfAS 

j 

O.01OA00 

n?  » 

0.0000*10 

0|  - 

7.  *160000 

n?  - 

*».  AAOOOO 

ir.  statssis 

(PM  » 

SHfAA 

SAIN 

Sf  r.HRNf 

-7767.0  -1 

7160.0 

-11*10.0 

v  t*uf  irf/srr.  i  * ao? . i 

f,  BA®  »  0.70A77ASF-06  PFO  HP'|* 

ANS 

Nr. 

f  IS10WI 

NC  FI  MOW) 

oust  vmir.iTY 

X  IAOOV 

if  mm 

1  Pf  A  I  T  1 

10. 

o. 

0. 1 1A 1A 1 If -OA 

0. 

I'. 

O.AWlosr-?? 

0.77*11A0n  -OA 

0.1ASAA77F-71 

?0. 

0.71*100161  -1  1 

0.0? 70OOHF-06 

0 .70  1  AA*76r  -  |  1 

77 . 

O.I*M?WPr-|| 

0;601 M667F-0A 

0.1 66  77  717  -  1 1 

?A. 

O.S7i|o/.or-|o 

0.1  7AA7 ’  If -06 

0.  MOO.ISOF-10 

76. 

,0.  aoh?  insr-nq 

0.7?  1 76  OOF  -  OA 

0.7AI 160 jr-oo 

7A. 

0.66007)  Al  -  OH 

0. 1  7 *76.161.  -  OA 

0.  I  777AA1  (j-OA 

10. 

0,  160in<}7(  -01 

0. 7  06  a 1 60r-0  7 

0. 16601 S1F-0H 

V. 

0.  1  A 1 A 1  0  Al  -Of, 

0.1771'  ’ Af -0 7 

0.7  7170S7F -OA 

1A. 

O.AS1*1A/0«  -06 

0.1«nS7AOf-0? 

0.17 167? IF -of 

1A. 

0.  1 1  nOAA  Tt  -OS 

0.01 AAAIor-OA 

0 , 1 S  7 AS Aflf -0 7 

1A. 

0.7AHnpi,f-0S 

0  ,  A  7-7  A  7  A  S  f  -  0  A 

O.16S*76O6f-07 

AO. 

0.S1O1  7MF-0S 

0.1*»J  177AF-OA 

O.IAOOIOM -07 

A6. 

0.  70AM  6 A 1  - 0 A 

0.71A0766F-00 

0.AAA01S6F-0A 

SO. 

0.  SA'JOAAlf  -OA 

0. 1 AS1 1 6  IF  - | 0 

0.1A67AA7P-0A 

ss. 

0.11  .’SHIP* -01 

O.I?!7?A«r--l  1 

0.700  7N01F-00 

60. 

0.IN1A1SAF-O1 

o.Aioo/,iv»r-i  i 

0.1 7AA6I IF- |0 

AS. 

0.7*»67'»‘1  ’f-io 

0 . 7 66 1  o/>  AF  - 1  A 

0.1001A04F-1J 

>0. 

0.  A1S0*»a’»  -01 

o. 7n 7o  \  >?r- 16 

0. A  7 1 7SAA7-I  1 

is. 

o.sasi  ii/.r-oi 

0.1*’  766  7  IF  -  J  7 

0. |A77?A6r-lA 

no. 

0. AM  no /S'  -01 

0.  1A6M  »AF-J*J 

0.  1.A0S  1  1  Af  -  16 

TABLE  20.  EXCEEDANCES  PER  HOUR  -  CONDITION  8  (33,000  FT. ) 


1  t*M  CONDITION 

§  VERTICAL  GUST 

ANALYSIS 

WING  CIA  STATION  •  0.11  STRUCTURAL  ELEMENT  110 

>  ALTITUDE  33*000  VT. 

TURBULENCE  PARAMETERS 

PI  -  0.012500 

P2  •  0.000090 

81  •  2.470000 

02  *10.400000 

1C  STRESSES  (RS1I 

SHEAR  SR  IN 

SEGMENT 

- 

2242.0  -12160.0 

-11810.0 

V  TRUE  IFT/SECt  >002.) 

G  BAR  •  0.4429345E-06  PER  HOUR 

KMS 

NC 

F  me  to 

NC  FI  SIGN) 

GUST  VELOCITY 

X  1800V 

(M/SEC) 

(PER  FT) 

10. 

0. 

0.  1594628E-04 

0. 

IS. 

0.4192 J95E-22 

0.2449926E-05 

0.1551958E-21 

20* 

0.2190014E-1 1 

0.1086651C-05 

0.1416544E-11 

22. 

0.191 2570E-1 t 

0.7369660E-06 

0.2035400E-IT 

24. 

0.57H040E-10 

0.48I6616E-06 

0.1986278E-10 

26. 

O.8OR2)05E-O9 

Oi  30)3  744E-06 

0. 154081 6E -09 

28. 

0.6600714F-08 

0. 1841420E-06 

0.I755209E-08 

10. 

0. 359)09 7E-07 

0.1077I24E-06 

0.5588814E-08 

12. 

0.143H104F-06 

0.6071 798E-07 

0.1260935E-07 

14. 

0.4539660F-06 

0.  32984 34F-07' 

0.2162100F-07 

16. 

0. 1 18964  7F-05 

0.1 726780F-07 

0.29666  72 C ”07 

18. 

0.26886  31E-05 

0.8711 743E-08 

0.3382171E-07 

40. 

0. 5  39  372  3E-05 

0.4235571F-08 

0.3299024E-07 

45. 

0.2086166E-04 

0.5938290E-09 

0.1 780938F-O7 

50. 

0. 549088  IE-04 

0.6607168F-10 

0.5239084E-08 

55. 

0.112  3830E-01 

0.58)462  7E-1 1 

0.9468880E-09 

60. 

0.1918156F-0) 

0.40B8984F-12 

0. 1 144410E -09 

65. 

0.2962882E-0) 

0.2274234F-11 

0.9730483E-11 

TO. 

0.4150862E-0) 

0. 1003859F-14 

0.601 7219E-12 

75. 

0.5451 336E-01 

0. 35166) 7F-16 

O.2760314E-11  j 

80. 

0.681 897HE-0) 

0.9776875E-18 

0.9627295E-15  1 
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TABLE  21.  EXCEEDANCES  PER  HOUR  -  CONDITION  9 


M0»  CUNOITIO*  6  VERTICAL  COST 

ANALYSIS 

NINO  If*  STATION  •  O.H  STRUCTURAL  ELEMENT  110 

TUR RULE MCE  RARAHETERS 

R|  •  0.00T000 

PI  •  o.oootoo 

01  •  J.0S0000 

02  •  6.1T0000 

1C  STRESSES 

IRSII 

SHEAR  SMN 

SEGMENT 

-166T.0  >12010.0 

>12260.0 

V  TRUE  (FT/SCCt  *002.) 

C  IAR  •  0.1665561E-07  RER  MCU« 

MS 

NC 

P  f Slew) 

nc  fisicmi 

CUST  VELOCITY 

i  iooov 

IfT/SECl 

IRER  Ml 

10. 

0. 

0.  D686S6E-06 

0. 

IS. 

C.2661067E-25 

0.2561715E-0S 

0.1 108 122E-26 

20. 

0.35062S6E-15 

0.6626S06E-06 

0.600S1TRE-IS 

22. 

0.6161S2TE-11 

0.  S06 1 866F-06 

0.S106501E-I) 

2A. 

0.  3263177E-U 

0.  16S6666E-06 

0.461027SE-I 1 

26. 

O.69E93STE-10 

0.16576011-06 

0. 16756I1E-10 

20. 

0. T660520C -06 

0.105R102E-06 

0.  12I660TE-06 

10. 

0.566602^-08 

0.S6M611E-0? 

0 .66566611-06 

12. 

0.2E12606E-07 

0.2667502E-0T 

0. 1065560E-08 

16. 

0. 1063STOr-06 

0.1272016E-0T 

0.16  51685E-08 

16. 

0.3262116E-06 

0. 57 11 1 2 l f -00 

0.2081 166E-08 

10. 

0. E12T062E-06 

0.26071 !!E“00 

0.2660668E-09 

60. 

0. 1063600F-05 

0.1016751E-00 

0. 2710561E-08 

61. 

0.06091665-05 

0.  6020050E-I 0 

0.11606TIS-00- 

SO. 

0.2T20T16E-06 

0.6011 1 1 7  E  —  1 1 

0.21T051TE-06 

ss. 

0.6267500E-06 

0.10111161-12 

0.2716A01E-10 

60. 

O.llTJItlE-01 

,0. 1167507  F- 11 

0.1666161E-U 

OS. 

0.16IOSOOE-01 

0..1205006E-IS 

0.6030I61E-11 

TO. 

0.20266766-01 

0.06601 7  1E-I 7 

0.205661 7E-14 

?s. 

0. 1R7O062E -01 

0.  1 12S.S50E-18' 

0. 6266666E- 16 

00. 

0.S016T12E-01 

0.1 1012 16E-20 

0.607571TE-18 

TABLE  22.  EXCEEDANCES  PER  HOUR  -  CONDITION  9  (38,000  FT. ) 

~  7208  CONOlVlON  9  VEAttCAL  GUST  ANALYSIS 

WING  ETA  STATION  •  0.11  STRUCTURAL  fLEMFNT  110 

ALTITUDE  38,000  FT. 


TUftSUlENCE  FAAAHETEAS 


El  •  0.008)00 
*2  •  0.000100 

81  •  2.00000 

82  •  8.500000 

to  STRESSES  1FSII 

SHE  A  A 

SKIN 

SEGMENT 

-1667.0 

-128)0.0 

-12260.0 

¥  TRUE  1FT/SECI  -802. 

) 

0  BAA  •  0.6021656E-07  EE  A  HOUR 

MS 

NC 

F  1 StGWI 

NC  FISIGWI 

OUST  VELOCITY 

*  1800V 

I F  T/ sec ) 

IEEA  FT  1 

10. 

0. 

0. 1200366E-06 

0. 

15. 

0.29936976-25 

0.2620095E-05 

0.1066225E-26 

20. 

0. 1589259E- l 5 

0.9157799E-06 

0.67665816-15 

22. 

0.6)61 52TE- l ) 

0.57501 83E-06 

0.52821606-13 

25. 

0.326)1 776-1 1 

0.1656012E-06 

0.162  77066-1 1 

26. 

0.69893576-10 

0.1986803E-06 

0.20032716-10 

28. 

0.79685206-09 

0.1091095E-06 

0.12523736-09 

)0. 

0.5669R25E-08 

0.5737982E-07 

0.6681636F-09 

)2. 

0.2812689E-07 

0.28867)06-07 

0. 11726196-08 

)♦. 

0. 1063578E-06 

0.1389)3^6-07 

0.21338366-08 

)6. 

0. 3262  31 9E  *06 

0.6396730E-08 

0.2995013E-08 

in. 

0.8327862E-06 

0.281  76  78E -08 

0.33882706-08 

60. 

0.1 861606E -05 

0. 11871 75E-08 

0.31965326-08 

65. 

0.89091 966-05 

0. 1 127036F-09 

0.16699776-08 

50. 

0.2728T15F-06 

0.61 10686E-1 1 

0.3195956E-09 

55. 

0.6267506E-06 

0.66266016-12 

0.39917726-10 

60. 

°  1 1 7)  31 3E-0) 

0. 1 829  7  32E-1 3 

0.31001766-11 

65. 

0.191 6566E-03 

0.57358506-15 

0.1S87673F-12 

TO. 

0. 28296 76E-0) 

0. 1 36)0206-16 

0.556962BF-16 

T5. 

0. 38  760626-03 

0.2655326E-18 

0. 1 376309F- 1 5 

80. 

0.5016712E-03 

0.3352615C-20 

0.26289226-17 

TABLE  23.  EXCEEDANCES  PER  HOUR  -  CONDITION  10 


7208  CONOITION  10  VERTICAL  GUST 

ANALYSIS 

** 

WING  ETA  STATION  •  0.31  STRUCTURAl  ELEMENT  110 

TURBULENCE  PARAMETERS 

PI  •  0.010500 

P2  •  0.000090 

81  •  2.950000 

82  •  9.850000 

1G  STRESSES  IASI  1 

SHEAR  SKIN 

SEGMENT 

-1600.0  -124*0.0 

-11R90.0 

V  TRUE  (FT/SEC)  -802.3 

C  BAR  «  0.942 185TE-07  PER  HOUR 

RMS  NC 

F  f SIGHT 

NC  FI  SIGH) 

CUSf  VELOCITY 

X  1800V 

t Ff/SECI  1  PER  FT! 

10.  0. 

0.1343433E-04 

0. 

15.  0.6825464E-24 

0.2293402E-05 

0. 2260463E-23 

20.  0.  190O235E  *'5 

0.9279008E-06 

0.2653404E-14 

22.  0. 2562953E-12 

0.6018662E-06 

0. 2227540E-12 

24.  0. 1034991E-10 

0.3746223E-06 

0. 5599060E-1 1 

26.  0.  1 839855E-09 

0. 223  7600E-06 

0.59449B8E-10 

28.  0. 1804632E-08 

0.  1282526E-06 

0. 3342257E-09 

30.  0.U38379E-07 

0. 7054160E-07 

0.U59624E-08 

32.  0. 5139145E-07 

0.  3723226E-07 

0. 2763091 E-06 

34.  0. 179214  7F.t06 

C.1885769E-07 

0.4880306E-08 

36.  0.5104520E-06 

0.9165430E-08 

0.6756048E-08 

38.  0. 1237806E-05 

0.4274765E-08 

0. 7640991E-08 

40.  0.2636902E-05 

C.  19 1 3228E-O0 

0.7285271E-08 

45.  0.  U45940E-04 

C.2140765E-09 

0.3542549E-08 

50.  0.3278019E-04 

?!.  1851253E-10 

0.8763187E-09 

55.  0. 71 351 82E-04 

0.123724RE-11 

0.1274814E-09 

60.  0.1289484E-03 

0.6390609E-13 

0. 1 I 099B9E- 10 

65.  0. 2044269E -03 

0.255106BF-14 

0. 7530B6BE-12 

70.  0.2S47557E-03 

0.7870392E-16 

0.3349991E-13 

75.  0.3961485E-03 

0.  11*763 71  t-l  7 

0.1073515E-14 

80.  0. 5049039E-03 

0. 345R024F-1 9 

0. 252 1 28  3E- 16 
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TABLE  24.  EXCEEDANCES  PER  HOUR  -  CONDITION  IT 

7200  CONDITION  It  V(MICU  CUSf  *N*l VSI S 
MINGE  TA  STATION  •  0.11  STRUCTURAL  flCNFNT  110 


lUARUtfNCC  FARAHFIERS 


Rt 

F2 

01 

02 

•  0.010900 

•  0.000090 

•  2.950000 

•  5.050000 

1C  SIRESSCS  IASI  1 

shear 

SKIN 

SECHfNI 

-1000.0 
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FIGURE  84.  SUMMARY  OF  DESIGN  VALUES  OF  EXCEEDANCE  RATIOS 
AND  FREQUENCY  OF  EXCEEDANCES  BASED 
ON  JOINT  PROBABILITY 


CONCLUSIONS 


General 


The  use  of  power  spectral  techniques  to  assess  the  strength  capability 
of  airplanes  subjected'  to  continuous  turbulence  appears  to  be  the  most 
rational  approach  to  the  -airplane  gust  design  problem.  The  statistical 
methods  which  were  developed  basically  in  the  communications  field  are 
well  founded  and  understood;  however,  their  application  to-  complex  air¬ 
craft  structure  in  dynamic  loads  analyses  has  only  been  pursued  in 
recent  years. 

There  are  two  rather  general  problems  associated  with  a  detailed  gust 
power  spectral  analyses.  The  first  is  the  "long  chain"  of  calculations 
that  must  be  conducted  to  arrive  at  the  solution,  and  the  second  is  the 
interpretation  of  the,  statistical  result- in  terms  that  are  .meaningful  to 
the  design  engineer,.  In  the  -foregoing  analysis,  the  first;  problem  was 
partially  overcome  by  proper  linking  of  the  individual  digital  computer 
programs  required  to  calculate  the  dynamic  responses  in  the.  flexible 
modes,  the  rms  loads  and  stivsses,  and  -finally  the  statistical  results 
into  a  chain  program.  This  allowed-  for  orderly  flow  of  data  by  the  use 
of  tape  storage  from  one  program  to- another..  This  enables  an  engineer 
to  analyze  many  conditions  in  a  relatively  -short  period  of  time.  The 
results  of  £his  study  will  serve  as  a  guide  to  the  solution  of  the 
second  problem. 

Dynamic  analyses  leading  up  to  an  assessment  of  the  strength  of  struc¬ 
tural  elements  or  components  is  complex  and  can  be  very  sensitive  to 
the  mathematical  description  of  the  airplane-  being  analyzed.  .'Even  when 
the  analysis  is  conducted  for  an  airplane  that  is  currently  in  service, 
a  vehicle  for  which  a  great  deal  of  information  is  available,  the  anal¬ 
ysis  is  still  difficult  to  accomplish.  Experienced  personnel  who  are 
capable  z±  exercising  good  judgment  on  estimates  of  sensitive  airplane 
parameters  are  needed,  to  use  this  approach  in  the  preliminary  design 
stage . 


Certainly,  power  spectral  gust  analyses  should  be  conducted  early  in  the 
design  stages  to  "raise-the-flag"  on  possible  problem  areas.  Further¬ 
more,  these  'analyses-  should  include  parametric  variations  of  sensitive 
parameters.,  which  need  better  definition.  The  analyses  may  become  more- 
complex  as  thf  design  nears  the  final  stages,  .but  the  analysis  should  be 
assessed,  ih  ?aqh  stage  ,by  comparing  reduced  .solutions  of  the  basic  equa¬ 
tions  with  separate  -Maneuver  and  flutter  solutions. 


183 


One  Parameter  Versus  Joint  Probability  Approaches 

If  the  structure  being,  analyzed  is  sensitive  primarily  to  an  axial  stress 
which  resu?-ts  predominantly  from  bending,  the  bending  moment  approach 
and  the  joint  probability  approach  should  lead  to  nearly  the  same  result. 
This  is  true,  because  the  other  load  quantities  contribute  only  slightly 
to  the  margin  of  safety  in  the  more  critical  structural  elements.  Since 
the  effort  involved  in  a  single  parameter  bending  moment  approach  is 
about  one-half  of  that  of  the  full  joint  probability  stress  analysis,  it 
is  advantageous  to  use  a  bending  moment  criteria  when  suitable. 

The  bending  moment  analysis  is  useful  in  the  .design  envelope  approach  to 
v  •'Ip  locate  the  critical  flight  condition  The  selection  of  critical 
conditions  for  vertical  analyses  in  this  program  required  20  separate 
analyses.  Ii  all  of  these  analyses  had  been  carried  through  the  joint 
probability  solution,  the  analyses  time  would  have  been  extended  several 
months.  Since  sill  of  the  load  quantities  tended  to  increase  as  the 
critical  condition  was  approached1,  the  bending  moment  served  as  .a  good 
measure  of  the  criticalhess  of  the  condition. 

The  success  of  using  the  bending  moment  to  locate  -the  critical  structure 
is  dependent  again  upon  whether  or  hot  the  structural  element  is  more 
likely  to  reach  limit  load  under  an  axial  stress  condition  or  a  combined 
tress  condition.  If  it  is  a  stress  condition  which  results  from  com¬ 
bined  loadings  then  the  bending  moment  alone  will  lead  to  erroneous 
results. 

The  joint  probability  method,  while  more  rational  in  considering  combined 
loading,  has  some  drawbacks  in  regards  to  its  use  in  a  design  sense.  One 
of  these  problems  is  in  obtaining  the  detailed  stress  coefficients, 
another  is  the  definition  of  the  strength,  envelope. 

Model  720B  Analysis  Results 

Limit -strength  values  of  orwT]^  and  of  hours  to  exceed  limit  design 
strength,  for  the  mode?  .720B,  are  summarized  in  Figures  85  and  86.  Re¬ 
sults  based  on  both  the  bending  moment  analyses  and  the  more  exact  joint 
probability  analyses  are  shown. 

Recommendations  for  Design 

Specifically,  the  following  recommendations  are  offered  concerning  the 
use  of  gust  power  spectral,  techniques  in  airplane  design: 

1.  Gust  power  spectral  analysis  techniques  appear  to  be  the  most 
rational  approach  to  the  gust  loads  design  problem.  However, 
they  should  be  used  in  such  a  manner  that  they  can  truly  as¬ 
sure  timely  design  information,  .and  hence  a  properly  designed 
airplane  from  the  standpoint  of  gust  loading.  This  can  be 


accomplished  by  a  continual  study  during  the  airplane  design- 
stages  using  the  best  .available  data  at  each  stage  to  eliminate 
the  possibility  of  unpleasant  surprises  after  the  airplane  is 
in  service. 

2.  initially,  the  aerodynamic,  structural,  and  stability  augmenta¬ 
tion,  data  used  in  the  analyses  will  be  rough;  therefore,  these 
initial  analyses  should  be  simple  and  should  include  parametric 
variations  on  the  more  significant  input  data. 

3.  Then,  the  more  sensitive  parameters  should  be  isolated  and 

*  studied  further  through  separate  test  and  analysis  programs. 

4.  Later  studies  will  be  Used  to  evaluate  the  limit  bending  moments, 
shears,  and  torsions  for  the  major  components.  Also,  the  phas¬ 
ing  or  correlation  between  these  loads  should,  be  obtained  and 
used  as  a  guide  in  the  design. 

5.  Joint  probability  techniques  to  properly  account  for  the  phas¬ 
ing  of  the  loads  or  stresses  should  be  employed  only  for  the 
.more  critical  components  of  the  airplane  when  the  margin  of 
safety  is  highly  dependent  on  combined  loading  and  the  statis¬ 
tical  correlation  between  the  loads  or  the  stresses  is  poor. 

6.  All  of  the  studies  should  be  subjected  to  a  set  of  checks  with 
separate  maneuver  and  flutter  analyses  or  test  results  in  order 
to  keep  the  studies  in  perspective  and  to  provide  confidence 

in  the  results-  as  the  analyses  become  more  complex  during  the 
later  design  stages. 

7.  Stability  augmentation  systems  and  control  dynamic  effects 
should  be  incorporated  in  the  early  studies  and  developed  a- 
long  with  the  basic  airframe  in  coupled  dynamic  loads  gust 
analyses. 
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AFPENDIX  A 
EQUATIONS  OF  MOTION 


General  Comments 

The  degree  of  success  in  evaluating  the  strength  capability  of  a  complex 
airplane  structure  subjected  to  flight  in  atmospheric  turbulence  relies 
almost  entirely  upon  the  ability  to  construct  a  mathematical  model  of 
the  physical  system.  An  airpiane  in  flight  is  considered  to  be  an 
elastic-dynamic  system  which  when  acted  upon  by  external  aerodynamic 
gust  velocities  undergoes  motions  necessary  to  maintain  an  energy  balance. 
The  analysis  procedure  generally  used  is  referred  to  as  the  Energy  Method 
and  is  applicable  to  a  system  where  the  total  change  in  kinetic  energy  is 
equal  to  the  summation  of  the  work  done  by  all  forces  acting  during  a  slight 
deflection  from  the  equilibrium  position-.  In  other  words  the  application 
of  the  energy  method  is  contingent  on  the  principle  of  virtual  work  which 
may  be  stated  as  follows: 

If  a  body  is  in  equilibrium  under  the  action  of  prescribed  external 
forces,  the  work  done  by  these  forces  in  a  small  additional  dis¬ 
placement.  compatible  with  the  constraints  is  equal  to  the  change  in 
strain  energy. 

The  above  statement  can  be  expressed  as  Eq.  (Al) 
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5w  +  5w,  «  5v 

e  m 


(Al) 


where 

5W  =  the  external  work  done 
*  e 

8 W.  =  the  internal  work  done 
in 

5V  =  the  change  in  strain  energy 

The  airplane  motions,  as  given  by  Eq  .yA2),  are  expressed  as  the  summation 
over  a.  set  of  independent  generalized  coordinates  which  are  defined  ex¬ 
plicitly  in  terms  of  space  and  time. 

n 

§i  =  Z  0ij  (A2) 

J=1 


A-2 


where 


<f>  . ,  =  the  jth  normalized,  mode  shape  which  is  a  function 
^  of  geometry  only. 

a  .  =  the  jth  generalized  coordinate  which  is  a  function 

^  of  time  only. 

Using  the  above  coordinate  definition  and  Lagrange's  specialized  form  of 
the  energy  equation,  the  equations  of  motion  for  the  physical  system  are. 
formed  using  Eq.  (A3). 
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where 


T  =  the  kinetic  energy  of  the  system. 
V  =  the  strain  energy  of  the  system. 


Qj  =  generalized  external  force  due  to  the  system's  response. 

Qj  =  generalized  external  force  due  to  the  gust  velocity. 

a  the-  jth  generalized  coordinate  displacement. 

=  the  jth  generalized  coordinate  velocity. 

For  a  linear  elastic  system  such  as  the  airplane  structure  the  elastic 
forces  are  given  by  d v/  and  the  kinetic  energy  is  not  dependent  upon 
the  displacement,  $j,  thus  dT/  5oi  =  0,  leaving  the  Lagrangian  equation 
as  Eq.  (A4)  below. 


VERTICAL  ANALYSES 


SL  |i 

i  'll 


Generalized  Inertia 

The  kinetic  energy,  T,  in  the  above  equation  is  obtained  by  multiplying  the 
individual  lumped  masses-  by  their  respective  velocities  squared  and  then 
summing  over  the  generalized  coordinates.  This  can  be  expressed  in  matrix 
form  as  Eq.  (A5). 
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where 


=  the  mass  and  inertia  properties  described  about  the 
principle  axis  at  the  lumped  mass  center  of  gravity. 

=  the  transformation  matrix  transforming  lumped  mass 
center  of  gravity  displacements  to  elastic  axis  dis¬ 
placements  . 

=  the  normalized  mode*  shapes. 


When  this  kinetic  energy  matrix,  T,  is  substituted  into  Lagrange's  equation 
and  the  appropriate  differentiations  carried  out  the  following  inertia 
matrix  is  generated: 
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V,  =  generalized  coordinate  acceleration 
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as  described  above. 


Generalized  Stiffness 

The  second  portion  of  equations  of  motion  by  the  Lagrangian  approach  are 
the  strain  energy  terms.  The  strain  energy  associated  with  a  deflection 
in  the  generalized  coordinates  could  be  generated  by  considering  the  elastic 
axis  as  a  series  of  springs  having  stiffnesses  of  k£.  Then  the  strain 
energy  would  be  given  by  Eq.  (AT ) * 

2v  ■  bllM'  M  M 


A-4 


and  when  substituted  in  Lagrange's  equation  would  result  as  Eq.  (A8) 
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where 


av 


=  kij|  ( "'ii]  Kj  | 


(A8) 


However,  it  is  more  convenient  to  use  the  natural  frequency  of  the  jth 
normal  mode  to  represent  the  generalized  stiffness  as  expressed  by  Eq.  (A9): 


A  tabulation  of  the  natural  frequencies  used  in  the  analysis  are  given 
in  Tables  6  and  7  in  the  Analysis-  Section. 

Generalized  response  Aerodynamics 

The  generalized  external  forces  which  result  from  the  motion  of  the  system, 
Qj,  are  the  aerodynamic  lift  and  moments.  These  force  expressions  are 
based  upon  a  modified  strip  theory  (°)'  which  incorporates  the  aerodynamic 
induction  between  the  individual  aerodynamic  panels.  The  transmittal  time 
associated  with  this  aerodynamic  induction  is  also  incorporated  as  a  phase 
lag  between  the  resultant  force  vector  and  the  originating  panel  motion. 

The  strip  theory  coefficients  which  give  the.  aerodynamic  forces  on  a  panel 
due  to  its  own  motion  are  used  as  they  were  originally  expressed  by 
Theodorsen  (9)  and  Kussner  (10)  with  two  exceptions.  These  exceptions  are 
the  two  dimensional  lift  and  moment  expressions,  2n ,  and  (l/2+a)  are  changed 
to  values  more  representative  of  the  aerodynamic  surface  being  examined. 

The  adjusted  values  are  obtained  by  reducing  the  oscillatory  lift  and  moment 
coefficients  to  their  steady  state  values  and  using  the  static  induction 
matrix  (?)  to  solve  for  the  lift  and  moment  coefficients  that  give  the 
proper  steady  state  lift  and  moment  distributions.  These  values  are  then 
used  with  the  dynamic  induction  matrix  and  oscillatory  lift  and  moment  co¬ 
efficients  to  give  the  generalized  external  aerodynamic  force  due  to  system 
response.  A  representative  equation  for  this  generalized  external  force 
is  given  below  in  matrix  form  as  Eq.  (A10): 
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where 
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2  4 

atmospheric  density  -  lb .sec.  per  in. 

reference  semi-chord  -  in. 

response  frequency  -  rad.  per  sec., 

span  of  the  ith  aerodynamic  panel  -  in. 

distance  from  panel  quarter  chord  to  panel  rotation 
axis  -  in. 

dynamic  induction  matrix  -  in. 

oscillatory  panel  lift,  coefficient  due  to  the  translation 
motion  'of  the  panel  rotation  axis . 

oscillatory  panel  lift  coefficient  due  to  angular  motion 
about  the  panel  rotation  axis. 

oscillatory  panel  moment  coefficient  due  to  translation 
motion  of  the  panel  rotation,  axis. 

oscillatory  panel  moment  coefficient  due  to  rotation 
motion  about  the  panel  rotation  axis. 


Selecting  a  frequence,  o>,  the  evaluation  of  the  above  generalized  external' 
force  expression  results  in  an  aerodynamic  matrix  with  complex  coefficients^ 
Since  this  expression  is  for  a  discrete  harmonic  frequency  these  complex 
coefficients  can  be  expressed  by  Eq.  (All). 


a+ib 


(All) 


where 
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i,  = 

2 


Applying  this  to  Q.  in  Eq.  (A10)  gives  a  set  of  aerodynamic  coefficients 
1  J,w&T 

for  the  generalized  coordinate  displacements  and  velocities  as  shown  in 
Eq.  (A12). 


(A12) 


The  aerodynamic  forces  acting  upon  the  body  panels  are  treated  separately. 
Quasi-steady  expressions  are  used  which  account  for  apparent  angles  of  at¬ 
tack  due  to  translation  velocities,  rotation  displacements  and  rotational 
velocities  of  each  body  panel.  The  expression  for  these  forces  are  given 
below  in  Eq.  (A13). 
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where 


\h 

\a 

Amh 

Ama 


quasi  static  lift  coefficient  due  to  translation 
quasi  static  lift  coefficient  due  to  rotation 
quasi  static  moment  coefficient  due  to  translation 
quasi  static  moment  coefficient  due  to  rotation 


These  are  calculated  for  each  frequency  and  converted  to  displacement  and 
velocity  coefficients  and  added  the  generalized  forces  calculated  for  the 
wing  and  horizontal  tail  to  give  total  generalized  external  force  due  to 
system  response . 


Generalized  Excitation  Aerodynamics 

The  remaining  external  force  arises  from  the  turbulence  velocities,  U , 

The  expressions  used  to  describe  this  generalized  force  are  similar  to 
those  used  to  describe  the  aerodynamic  forces  due  to  response.  The  only 
essential  difference  is  that  these  forces  arise  from  gust  angles  of  attack. 
When  employing  the  power  spectral  density  approach,  it  is  necessary  to  ob¬ 
tain  the  response  of  the  airplane  to  discrete  gust  frequencies.  These 
gust  angle  inputs  are  expressed  by  Eq.  (Al4). 
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gust  velocity,  usually  taken  as  1  ft.  per  sec. 
true  airplane  velocity,  ft.  per  sec. 
gust  excitation  frequency,  rad.  per  sec. 


(Al4) 
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As  the  airplane  penetrates  the  gust,  the  excitation  force  is  felt  progres¬ 
sively  along  the  direction  of  flight.  The  gradual  penetration  effect  is 
represented  as  a  phase  lag  by  Eq.  (A15). 


where 
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the  distance  from  the  leading  edge  of  the  foremost  panel 
to  the  mid-chord  of  panel  i. 


Substituting,  the  expression  for  $  into  Eq.  (A15)  gives  Eq.  (Al6) . 
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The  end  expression  for  the  generalized  external  force  due  to  a  sinusoidal 
gust  velocity  is  given  by  Eq.  (A17). 


where 
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modified  two  dimensional  lift  coefficient  due  to  gust 
angle  of  attack. 


Structural  Damping 

The  final  item  that  may  be  included  and  is  not  shown  in  Lagrange's  equa¬ 
tion  is  the  structural  damping  force.  This  force  is  assumed  to  be  pro¬ 
portional  to  the  generalized  stiffness  but  in  phase  with  the  generalized 
coordinate  velocity.  An  expression  for  this  is  shown  below: 


(Al8) 


structural  damping 


1  gj  kj  % 


where 


damping  coefficient  approximately  equal  to  twice  the  damp¬ 
ing  ratio  c/c  . 

phasing  operator  which  shifts  the  force  vector  90  degrees 
ahead  of  the  displacement 

generalized  stiffness  for  the  jth  mode. 


The  expression  used  in  the  equations  of  motion  is  given  by  the  following 
form:  „  2 

1 kj  v  '  — sr-2 —  fa.  (A19) 


The  equations  of  motion  for  the  mathematical  model  of  the  airplane  system 
are  formed  by  collecting  all  of  the  above  expressions  into  a  matrix  equa¬ 
tion  as  follows: 
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or  simplified 
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Substituting  for  5.  and  5.  under  the  assumption  of  harmonic  motion  as  shown 
in,Eq.  (All)  we  obtain 
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Solution  of  Equations  of  Motion 


The  final  step  is  the  solving  of  the  equations  of  motion  to>  obtain  the  re¬ 
sponse  of  the  generalized  coordinate, .  This  solution,  obtained  for 


each  discrete  frequency,  6)  as  shown  by  Eq,  (A23),  results  in  the  vector 
magnitude  and  phase  angle  relative  to  the  gust  input  of  each  of  the  gen¬ 
eralized  coordinates. 
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Representative  frequency  response  curves  are  shovn  in  Fig.  2 J  of  the 
Results  section. 


LATERAL  ANALYSES 

Areas  of  Similarity  between  Vertical  and  Lateral  Analyses 

The  equations  of  motion  for  the  lateral  analyses  are  also  formed  using  the 
Lagrangian  approach.,  The  formulation  of  the  generalized  inertia  and. 
stiffness  matrices  for  the  lateral  equations  of  motion  follow  the  same 
procedure  as  outlined  for  the  vertical  analyses.  The  generalized  coordi¬ 
nates  used  in  these  lateral  equations  are  the  antisymmetric  body  and 
vertical  tail  cantilevered  modes  as  listed  in  Table  8  of  the  Analyses 
Section. 


Generalized  Aerodynamics 

The  aerodynamic  response  and  excitation  forces,  however)  axe  based  on 
quasi-steady  expressions  which  are  convoluted  with  the  Wagner  and  Kussner 
lift  growth  functions.  The  Wagner  and  Kussner  lift  growth  effects  were 
represented  by  the  step  response  functions  (recommended  on  Page 
Ref.  18)  for  an  aspect  ratio  of  three.  The  frequency  domain  analysis  re¬ 
quires  that  lift  growth  effects  be  represented  by  the  impulse  function 
which  is  obtained  by  writing  the  Laplace  transform  of  the  first  time  de¬ 
rivative  of  the  step  function. 


Yaw  Damper  Simulation 

The  yaw  damper  is  designed  to  reduce  'the  yaw-  response  to  disturbances  and, 
in  doing  so,  increases  the  Dutch  roll  stability.  The  yaw  damper  system  is 
best  defined  as  a  electromechanical  and  composed  of  a  rate  gyro  which  pro¬ 
vides  an  output  signal  proportional  to  yaw  rate.  The  signal  from  the  rate 
gyro  is  applied  to  a  bandpass  filter  which  removes  the  low  frequency  com¬ 
ponents  resulting  from  steady  state  turns  and  the  higher  frequency  compo¬ 
nents  associated  with  the  structural  vibrations.  After  filtering,  the 


signal  is  amplified  and  applied  to  a  servo-motor  which  results  in  a 
hydraulically  actuated  rudder  deflection  to  oppose  the  yawing  motion. 

The  method  of  simulating  the  yaw  damper  is  illustrated  by  the  simplified 
block  diagram  below: 


where : 

=  yaw  rate,  deg/ sec. 

eQ  <b  input  to- summation  point  from  filter. 
5gv  =  servo  output-drum  rotation,  deg. 

*  rudder  rotation,  deg. 


The  transfer  function  of  the  servo  amplifier  and  motor  combination  is  ap¬ 
proximated  by 


The  response  time  of  the  servo  is  assumed  to  be  negligibly  small  and  the 
expression  for  e0  can  be  expressed  by  the  following  equation. 


A-ll 


Solution  of  Equations  of  Motion 


The  input  format  for  the  digital  computer  program  requires  that  coeffi¬ 
cient  matrices  he  formed  as  indicated. 


hi  vKUrhU* -  hi*'* 


where : 


* 


generalized  stiffness  matrix. 

noh-circulatory  response  aerodynamic  damping  matrix, 
generalized  inertia  matrix. 

circulatory  response  aerodynamic  stiffness  matrix, 
circulatory  response  aerodynamic  damping  matrix, 
generalized  excitation  aerodynamic  matrix, 
instantaneous  gust  angle. 

indicates  the  application  of  Duhamel's  integral  to 
include  Wagner  or  Kussner  effects. 


The  generalized  force  coefficients  associated  with  rudder  deflection  due 
to  the  yaw  damper  are  formed  into  a  set  of  matrices  defined  as  j  SjT(s) I 
in  the  digital  computation  operation  format.  The  coefficient  matrices* 
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for  the  equations  of  motion  are  then  modified  from  the  above  form  to  in¬ 
clude  the  yaw  damper  effects  as  shown  by  the  following  expression: 
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where 


M. 


M. 


•the  new  coefficient  matrices  of  <^}  ^  and 
the  original  coefficient  matrices  of  fji  and 


The  digital  program  used  for  the  lateral  solution  solves  for  the  responses 
by  using  the  Laplace  transformed  equations  of  motion  with  the  Laplace 
transform  variable  set  equal  to  io.  The  solution  then  proceeds  as  follows: 
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^  ■  Laplace  transform  of  1st  time  derivative  of  Wagner 
f  =  Laplace  transform  of  1st  time  derivative  of  Kussner 


This  procedure  yields  a  set  of  complex  generalized  responses  for  the  fre¬ 
quency  under  consideration.  The  complex  generalized  responses  are  then 
used  to  compute  the  load  responses.  A  set  of  representation  solutions  for 
for  lateral  Conditions  10  and  10YD  are  given  in  Pig.  28  of  the  Analyses 
Section. 
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Crossings  of  an  Arbitrary  Strength  Envelope 
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STATISTICAL  ANALYSIS  FOR  COMBINED  RANDOM  STRESSES' 


Design  Envelope  Analysis .Stress  Relationships 

The-  cumulative  probability  function,  P(MS  <  0),  is  derived  from  the 
strength  enveippe  for  the  structural  element  in  question,  the  Gaussian 
joint  probability  density  distribution  for  the  incremental  axial  and 
shear  stresses  on  the  structural  element;  and  the  one-factor  level  flight 
stresses.  Consider  the  diagram  in  Fig.  Bl  where  'the  joint  probability- 
density  function  for  axial  and  shear  stress  is  shown  in  relation  to  the 
strength  envelope  defined  by  the  appropriate  interaction  curves.  The 
probability,  P(MS  >  0,  aw),  that  the  combined  axial  and  shear  stresses 
wiil  not  exceed  the  limit  strength  for  a  given  root-mean-square  gust 
intensity  is  equal  to  the  volume  of  the  -probability  function  within  the 
region,  R,  bounded  by  the  appropriate  limit  strength  envelope. 


P(MS  >  0)  =  ^p(f,£)ufd£ 


(Bl) 


Of  course,  the  probability  that  the  margin  will  be  less  than  zero,  that 
is,  that  the  margin  will  be  negative  is, 


P(MS  <  0 ,aw)  =  1  -  P(MS  >  (  ,aw) 


(B2) 


The  probability  function,  P(MS  <  0,aw),  can  be  obtained  by  several  inte¬ 
grations  for  various  aw  and  plotted  as  illustrated  in  Fig.  B2. 


L. 


The  joint  Gaussian  probability  density  function  illustrated  in  Fig.  Bl 
can  be  expressed  as  follows: 
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where 


1  T(f-f)2  {(-l)2 

(1-P2)  at  of 
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The  parameter y  P ,  is  the  correlation  coefficient  expressing  the  statis¬ 
tical  correlation  between  the  axial  and  shear  stresses.  It  can  be  de¬ 
fined  from  the  cross-power  spectrum  between  the  axial  and  shear  stresses 
as  follows: 


PofOg=  f^((o)d<o  =  ave  jf(t)f(t) 


Eq.  (b4)  is  derived  as  follows:  Assume  that  hf(t)  and  h£(t)  are  im¬ 
pulse  response  functions  for  the  axial  and  shear  stresses.  Then,  the 
time  histories  of  axial  and  shear  stresses  for  an  arbitrary  gust  loading, 
w(.t).,  can  be  determined  by  convolution. 


f(t)  =  jfhf( 
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£  (t-r)  »  f  h^(r2)w(t-r-r2)dr2 
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By  definition  the  cross-correlation  function  for  the  axial  and  shear 
stresses  is, 


Rf  (r)  =  lim  i  J  f(t)f(t-r)dt. 
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If  we  assume  that  the  gust  loading  has  been  acting  prior  to  time  zero, 
the  lower  limits  on  the  convolution  integrals  can  be  extended  to  minus 
infinity;  Also,  the,  order  of  integration  can-  be  changed. 


Rf*(0  =  f  /hf(r i)h>(r2)  lim  /Sw(t-r1)w(trr-r2)dt  d^dfg  (B7 
C  -J00  -x  * 
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The  cross-power  spectral  density  function  is, 


=  £  /Rfe(r)e-ia,rdr 
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If  we  set  r.-r^+fg  *  r^  anc*  rewrite  the  interior  integral  we  obtain, 
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Rearranging  again,  we  obtain 
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However,  we  note  that 
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where, 


0W(*>) 

Hf(«) 


gust  power  spectral  density  function 


complex  frequency  response  function  for  axial  stress 
in  the  structural  element 


h£(*>) 


complex  conjugate  frequency  response  function  for 
shear  stress  in  the  element: 


Therefore, 


0  =0w(")Hf(«)H|(«) 


(b: 


Obviously,  the  cross-power  spectral  density  function,  is  complex,  that  is 


0  =  Cf^6>)  +  iqf^(w) 


(Bll.) 


where, 


Cff(«) 


real  part,  cospectrum 

imaginary  part,  quadrature  spectrum 


Also  we  note  that, 


=  c^f(6>)  =  C^-Ol)  (B12a) 

qf^)  =  »  -qf^(-w)  (B12b) 

Therefore, 
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and 


oo 


|  /0w(&>)Hf(w)H>(w)dO> 


-oo 


r  oo 


/<^w(«»)|Hf(«)|  d6>  /^w(6>)|h.(0))|  dto> 


(B13) 


(  Applied  shear  stress 

MS  Margin  of  safety 

Setting  the  margin  of  safety  equal  to  zero  and  solving  for  the  normal 
stress  in  terms  of  the  shearing  stress  one  obtains 
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The  resulting  envelope  is  shown  in  Fig.  B3« 

Flight  -  Prof ile  Analysis  Stress  Relationships 

The  number  of  times  per  hour  that  zero  or  negative  strength  margin  will 
occur  on  a  structural  element  for  a  given  flight  condition  can  be  ex¬ 
pressed  as  follows: 
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where, 

f(aw)  the  probability  density  distribution  of  rms  gust 
''  velocity  for  the  appropriate  altitude  is  given  by 
Eq.  (2  ),  and  is  repeated  below 


(B22) 


P2,  b^,  and  bg  are  those  given  in 


The  constants  Pl' 
Figs.  5  and  6. 


Nc(°w)  Number  of  crossings  of  the  strength  envelope  per  second 
for  two  correlated  random  stress  time  histories.  The 
derivation  of  Nc(<rw).  is  given  in  the  following 
discussion. 

Crossings,  of  an  Arbitrary  Strength  Envelope 
Consider  a  probability  density  function 

'P ( x > a > y>/3) axdadyd/3  =  Prob  [ x<£<x+dx , a<£< a  +da,  y<f<y+dy, 

0<f</9+d0] 


where, 


f  axial  stress,  psi 

/3=  f  time  rate  of  change  of  eixial  stress,;  psi  per  sec. 
£  shear  stress,  psi 

a  «*  g  time  rate  of  change  of  shear  stress,  psi  per  sec. 


For  a  given  period  of  time,  the  above  expression  rep£esents  the  fraction 
of  time  per  unit  time  that  the  stress  vector  z’(t)  -  x(t)i  +  y(t)j  spends 
in  the  interval  (x,x+dx;  y,y+dy)  with  velocities  in  the  intervals 

(a,a+da;  /3,/3+d/3) . 

The  expected  number  of  times  per  unit  time  that  the  stress  vector,  'z(t), 
passes  through  the  interaction  boundary  interval  (Zf,(x,y),  z^(x,y+dij)  is 
equal  to  the  time  spent  in  the  interval  divided  by  the  time  required  to 
cross  the  interval. 

The  time  required  to  cross  the  interval,  dij,  is  equal  to  d  17  divided  by 
the  velocity  component  of  “  normal  to  C,  the  interaction  boundary  curve 
indicated  in  Fir.  b4. 

A  unit  vector  normal  to  C  is, 


grad  C 
|  grad  C 


(B24) 


If  C  is  a  function  of  x  and  y,  such  as 


where  this  function  represents  ah  interaction  curve  for  shear  and  tension^ 

x2 

C .  -  y  -  Ft  +  f~  -  .0 
t 

The  gradient  of  C,  denoted,  grad  C,  is 


Then,  the  rate  at  which  the  interval  d»?  is  traversed  is  equal  to  the  dot 
product  of  the  unit  normal  vector  with  the  velocity  vector. 


+ 
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T/ 


,/3?t 

v/ilx2  +  F| 


If  the  time  required  to  cross  the  interval  dij  is  f,  then. 


dii 

r  =  - - — - - 

grad  C  ,  i 

grad  C  l 


(r  >0) 


(B26) 


The  expected  number  of  times  per  unit  time  that  the  curve  C  is  crossed 
by  T  with  velocities  a  and  /8  is  equal  to  the  time  spent  in  the  interval 
d i\  divided  by  the  time  required  to  cross  the  interval. 


[E(x,a,y,/3)dad/3]dij 

I  grad  C.  ±.\ 

•  • "  * . * . •  z 

I  | grad  C|  1 


dS 


(B27) 


The  total  number  of  passages  both  into  the  interval  d»j  from  the  region  R 
and  into  the  interval  toward  the  region  R  is  obtained  *by  integrating  for 
all  possible  velocities  a  and/3.  Therefore,  the  total  number  of  pass¬ 
ages,  Nc  is, 


n*  -  /  /!/[— 

p  =-oo  a--  oo  i  1 


dsf  dad/3  (B28) 


Now,  if  the  process  is  Gaussian,  the  probability  density  can  be  written 
as  follows: 


where  |  M |  is  the  determinate  of  the  time  averages  of  the  products  and 
cross  products  of  x-^XgjX^,  *  *  * ,  x^,  and  M-j  is  the  co-factor  of  the  ij 
element  of  |m|. 

If  we  let 

Xjl  =  x(t) 

x2  ■  ®(t)  =  x(t) 

x3  -  'y(t) 

xif  -£(t)  =  y(t) 


then  |m| 


can  be  expressed  as, 
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T?' 

x(t)‘ 

[x(t)  o(t)  y(t)  /3(t)J  dt 

1] 

O(t) 

' 

y(t). 

J. 

T 

(B30a) 


or 


f 

r  2 

X 

xa 

xy 
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y£ 

-0X 

/3a 

fly 

/32J 

dt 


(B30b) 


The  elements  of  |M|  are  determined  b?  carrying  out  the  integration  over¬ 
time  to  obtain  the  following  for  the  (i,i)  elements. 
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r  li«i~  J  x2(t)dt 
x 


<*v 


X  *3 

=  a2  =  lim  ^  /  a2(t)dt  =  f  u^(<u)dw 


T-**-oc 


X  w 

=  a2  =  lim  ~r  y*  S^UJdt  =  ^  £y(6>)d« 


T-^oc 


°I  -  „lim  4  /  p2(t)< 

M  T  -*■  oo 


w 

-  /4wi“ 


For  the  (l,2)  and  (2,l)  elements  we  obtain, 


a 

lim  5?  /  x(t)i(t)dt  = 


T-*-oo 


I |  !  This  is  determined  by  integrating  by  parts  and  by  taking  the  limit 


T  |  T  T 

lim  ~  j  x(t)x(t)dt  =  lim  ^ !  [x(t)x(t)J  -  f  x(t)x(t)dt 

'•^x  _T  T-^-oo  |  jj>  _<p 


Then, 


lim  i  /  x(t)x(t )dt  =  lim  ^  Tx(t)  -  x(-T)l  =  0 

T«^oo  _tp 


This  integral  is  identically  equal  to  zero,  since-  the  respective  values 
of  x(T)  and  x(-T)  are  always  finite  and  they  are  divided  by  T  which  is 
arbitrarily  large. 

Therefore,  the  (l>2)  and'  (2,l)  elements  of  |m|  are  zero. 

The  (l,3)  and  (3>l)  elements  are  by  definition  of  the  cross  correlation 
between  x(t)  and  y(t)  and  are 


f  X 

xy  =  Poa  =  lim  ^  J  x(t)y(t)dt  =  | /*^xy^dw  (B32). 
y  T«^oo  1 


where  p  is  the  correlation  coefficient  which  varies  in  the  interval 

(1,-1). 

All  other  elements  are  more  obscure  will  be  assumed  equal  to  zero. 


Therefore, 


|M|  = 


°x  0  P*x<fy  0 

0  o'i  0  0 


P^xfy  0  a2  0 

J 

0  0  0 

P 

The  frequency  density  is  now  expressed  as  follows: 


(B33) 


P(x,a,y,/3)  = 


(2ff)2  Va^(l-^)1/2  I2(l-p2)|a2  axay  a2 


-1  jx2  2pxy  y2 


♦4 


(B34) 


If  we  assume  that  x  and  y  shown  in  Eq.  (B34)  are  incremental  or  time 
variant  values  varying  from  the  mean  or  one-factor  flight  values,  then, 
x  is  replaced  by  (x-x0)  and  y  by  (y-y0). 
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p(x,«,y,/3)  = 


1 _ l  j  .  -l  Xx-*cf 

(2f)2  <7x0’0ffy^(l-p2)1/5  (Sd-P2)! 


'0/3  (2n): 


i*“T»pr3F0TOT 


expU!4+4l!/1^44^ 


(x-x0)(y-y0)  (y-y0) 


£44#  411 
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Substituting  Eq.  (B35)  into  Eq.  (B27),  we  obtain 


(B35) 


-1  [(x-x0)2  o  (x-x0)(y-y0)  ( y-y0 f  ,a0  ... 

-i'2— [— — *  <B36> 


The  number  of  crossings  per  unit  time  is  obtained  from  Eqs.  (B28)  and 
(B35). 


tap*  fit 


1 _ 1  r  f  rl  grad  C  i 
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(B37) 

+fi~]]dSdad/3 


If  we  assume  zero  mean  values  for  x(t)  and  y(t),  we  can  easily  obtain  the 
familiar  expression  for  the  number  of  times  per  unit  time  that 
*2r(t)  =  x(t)i  +  y(t)j  crosses  a  constant  level,  y  =  a. 


**»r  - 


The  general  equation  for  the  number  of  crossings  per  unit  time  is  ob¬ 
tained  in  Eq.  (B37).  If  we  assume  that 


C  =  y -  aj 


9i  • 0i 


f-1 


and  the  grad  C  «  ,0i  +  lj.  Then,  (grad  cj 
component  of  ai  +  ft  j  normal  to..  C  is, 


=  1.  The  veloi 


/3=  (Oi  '41/3,))  •  (ai  +  /3a) 


Then, 


grad  C 


| grad  C 
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If  Ea.  (B37)  is  rewritten  for  this  example,  we  obtain, 
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>  clad/3 


00  I 

-1  r: 

1  2 it  p.1/8  exp 

»x<VCi-o2)  1 

’X  Vf  < 


We  note  that  integral  over  x  is  identically  equal  to 
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/■ 


x=-» 


v/2W, 


ity  vector 


(B38) 

|  dx 


(B39a) 


and  the  double  integral  can  be  rewritten  in  simpler  form  after  the  follow 
ing  integrations  are  made: 


/  e  da  =  v  2»r  crfl 


since. 


a  =-oo 


(B39b) 
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(B39c) 
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The  integral  over  /3  can  be  written  as, 


/  1*1 

fi=~oo 


-f?  /2<>p 


-  2  I  Be 


The  double  integral  over  a  and  /3  is  equal  to 

//■ 

«  |9 

Substituting  Eqs.  (B39a)  and  (B39b)  into  Eq.  (B38)  gives, 


(B39d) 


(B39e) 


Nr  -  {2\fzn  on o%)  ~== 

C  2”oaop  v  a  /3  v/15^( 


2,  1  -y2/2ay 


Reducing  Eq.  (B39f)  gives, 


-y^/Sc 


(B39f) 


(bUo) 


This  is  the  well  '-nown  expression  for  the  total  number  of  crossings  of  a 
level  y  =  a,  where  a  is  arbitrary  in  magnitude  and  in  sign.  C11) 


The  preceding  derivation  treated  a  continuous  envelope  C;  however,  the 
algebra  is  simpler  if  the  envelope  is  defined  by  a  series  of  straight 
line  segments.  In  the  following  development,  an  expression  for  the  value 
of  Nq  is  determined  for  one  of  these  line  segments.  Through  a  cycling 
process  the  Nq  for  any  envelope  can  be  determined. 

Consider  the  number  of  crossings  per  unit  time  of  an  arbitrary  straight 
line  segment  from  (x^y^)  to  (xqjYq)'  The  equation  of  the  line  through 
the  point  is, 


y 


(x-Xj^)  +  yx 


(B4la) 


or 


y  =  n^x-x^  +  y1 


where  m  is  the  slope. 
Therefore, 


C  -  y  -  n^x-Xjj  -  yx  -  0 


^2  =  -m;  -  lj  grad  C  =  -mi  +  j  (b42) 

Ox  O  y 


The  absolute  value  of  the  grad  C  is, 


grad  C 


l| 
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-ma+B 

>/m^  +  1 


Assuming  xo  =  y0  =  we  can  rdwrd-te  Eq.  (B37)  as, 


"c  ’  (*->2  «Av^)] 


r  r  f  1  c:.J  -1  fx2  _  +  yfl ) 

f  f  I  \fiBi  (2d-^)k  »2J'! 


dS  (b44) 
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Ea.  (b44)  can  be  rearranged  and  applied  to  the  k-th  segment  of  the 
strength  envelope,: 
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(b45) 
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xl>yl 

Eq.  (B45)  shows  that  Wc  is  the  product  of  a  slope-velocity  function, 
which  is  constant  for  a'^given  line  segment,  and.  a  line  integral  over  the 
Joint  probability  function,  P(x,y).  While  it  is  possible  to  express 


both  functions  in  series  form,  it  has  been  found  to  be  ^advantageous  to 
evaluate  the  slope-velocity  term  by  an  infinite  series  and  to  evaluate 
the  second  function  directly. 


The  suggested,  method  of  solution  is  to  numerically  evaluate  the  crossing 
density,  dNq  /dS,  at  a  sufficient  number  of  points  along  the  various  line 
segments  around  the  strength,  envelope. 
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The  crossings1  of  the  k-tb  segment  of  the  strength  envelope  can  be  de¬ 
termined  as  follows: 


Eq.  (b47)  can  be  properly  evaluated  if  the  series  converges  and  if  the 
second  term,  the  line  integral,  is  carefully  evaluated  ever  each  line 
segment.  The  convergence  of  the  series  in  less  than  30  terms  can  be  as¬ 
sured  if  (l)  the  ratio  y°)3  is  normalized  to  unity,  and  if  the  line 
segment  slopes  are  restricted  to  the  interval  (-1,1). 

The  ratio  °a/aB  can  be  made  equal  to  unity  by  stretching  the  strength 
envelope  and  all  x-dependent  parameters  by  the  ratio  V<jg.  The 
stretching,  of  course,  will  change  the  slopes  of  the  line^ segments.  Nc^ 
can  be  evaluated  directly  after  transformation  for  all  line  segments 
with  slopes  in  the  interval  (-1,1).  Segments  with  slopes  outside  the 
interval  (-1,1)  can  be  handled  by  exchanging  the  coordinates  in  Eq.  (b47). 
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Combined  Stresses,  for  Vertical  and.  Lateral  Gusts 

Calculations  to  estimate  airplane  loads  and  stresses  resulting  from 
flight  in  continuous  atmospheric  turbulence  are  nearly  always,  carried 
out  under  the  assumption  that  the  vertical  and  lateral  motions  of  the 
airplane  are  not  coupled.  As  a  result,  two  separate  analyses  are 
usually  conducted;  one  for  the  vertical-  component  of  the  turbulent  gust 
velocity  field  and  a  separate  analysis  for  the  lateral.  However,  cer¬ 
tain  parts  of  the  airplane,  -particularly  the  fuselage  and  empennage,  are 
stressed  significantly  by  the  simultaneous  action  of  both  the  vertical 
and  lateral  gust  components.  It  is  shown  in  the  following  development, 
that  stresses  derived  separately  from  vertical  and  lateral  power  special 
gust  analyses  can  be  combined  in  a  simple  manner  if  the  turbulent  gust 
velocity  field  is  isotropic. 

Consider  two  stress  components  in  a  structural  element  in  the  fuselage. 
Let  fv(t)  be  the  axial  stress  time-history  that  results  from  the  ver¬ 
tical  component  of  turbulence,  and  let  (t)  be  the  shear  stress  re¬ 
sulting  from  the'  lateral  component  of  turbulence.  The  cross  correlation 
between  these  two  stress  components  can  be  expressed  as  follows: 


T 

Rf  c  ( 0  *  lim  4  f  fv(t)&(t-r)dt  (B48) 

v$L 


Expressions  for  the  time-histories  of  fv(t)  and^^(t)  can  be  written  in 
terms  of  the  stress  impulse  response  functions  hfv(  r)  and  h^(  f)  and  the 
gust  time-histories  using  summation  integrals. 


fv(t)  =  J  hf^(rL)w(t-r1)dr1  (B49a) 

-  oo 
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Eq.  (-B48)  can  now  be  rewritten  as  follows: 
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-T  L' 


=  „lim^  /  /  hf 

_  “oo 
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jf  h^(r2)v(t-r-r2)dr2 

-  bo 


dt 


Rearranging  the  order  of  integration,  we  obtain. 
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(B50b) 


It  will  be  noted  that  the  term  in  brackets  in  Eq.  (B50b)  is  the  cross 
correlation  function  for  the  vertical  and  lateral  gust  velocity- 
components. 


W-v-'s) 


T 

lim  /  w(t-r1)v(t-r.r2)dt 
^"♦■00  -T 
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The  cross-power  spectral  density  function  for  the-  two  stress  components 
is, 


*;  fa* 

•00 
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where, 


=  r  -  r,  +  r , 


(B53) 
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After  substituting  Eqs.  (B53)  and  (B5l)'  into  (B50b)  and  rewriting  Eq. 
(B52),  we  have, 
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(B54). 


Rearranging  the  order  of ' integration  in  Eq.  (B5*0  gives. 


where. 


H,  (to) 


6  (a>) 

vwv'  7 


complex  frequency  response  function  for  axial  stress  on 
the  element  resulting  from  the  vertical  component  of 
atmospheric  turbulence 

complex  conjugate  frequency  response  function 
for  shear  stress  on  the  element  resulting  from  the 
lateral  component  of  turbulence. 

cross-power  spectral  density  function  for  the  vertical 
and  lateral  gust  velocity  components. 


The  integral  of  the  crossrpower  spectrum  in  Eq.  (B55)  gives 


pf  C  Of  <7> 
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where  pr  p.  is  the  correlation  coefficient. 
v5L 


Eq.  (B56)  provides  the  same  result  as  the  following  integral: 
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Pf  g  ~  li,n  ^[T  /  (B57) 
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The  axial  and  sheaf  stresses  -at  a  -.point  in  the  airplane  structure  re¬ 
sulting  from  the  vertical  component  and  the  lateral  component  of  atmos 
pheric  turbulence  are  as  follows: 


f(t)  =  fv(t)  +  fL(t) 

f(t)  =ev(t)  +fL(t) 


’(B58) 


The  matrix  of  mean  squared  stresses  used  to  define  the  joint  probability 
distribution  is. 
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(B59) 
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If  it  is  assumed  that  the  turbulence  as  isotropic,  the  time  averaged 
product  of  two  .gust  velocities  w(x^,t)  and  v(x2,t)  at  two  points  x^  and 
Xg  along  the  airplane  fuselage  is, 
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It  would  be  expected  that  the  above  space  correlation  would  be  sensitive 
to  the  absolute  value  of  t  minus  Xg.  However,  in  the  same  manner  as 
von  Karman  and  Howarth'1-5',  it  can  be  shown  that  this  correlation  is 
always  identically  zero  for  isotropic  turbulence  for  any  x^  and  Xg. 

If  we  rotate  the  coordinate  axis  system  about  either  the  y  or  z-axis  by 
].80°,  and  denote  the  transformed  velocities  and  v',  respectively,  we 
observe  that 


wv  a  -wv"  (B6la) 


But  by  the  isotropic  property. 
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Therefore, 
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Since, 
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then, 
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Therefore,  the  vertical  and  lateral  components  of  isotropic  turbulence 
are  not  correlated,  and  the  cross  power  spectral  density  function  for 
the  vertical  arid  lateral  components  of  gust  velocity  is  identically  zero 
for  ail  frequencies.  Now,  by  Eq.  (B5 6).  we  see  that  ail  of  the  correla¬ 
tion  coefficients  relating  stress  components  arising  separately  from 
the  vertical  and  lateral  components  of  isotropic  turbulence  are  also 
equal  to  zero. 

The  matrix  of.  mean  squared  stresses  used  to  define  the  joint  probability 
distribution  is  obtained  by  integrating  the  remaining  terms  in  Eq.  (B59)* 
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Now,  the  total  mean  squared  stresses  from  the  vertical  and  lateral  analy¬ 
sis  results  are, 


(B64) 


and  the  correlation  coefficient  relating  the  total  axial  and  total  shear 
stresses  is. 
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The  power  spectral  density  function  for  the  total  axial  stress  can  be 
obtained  from  the  total  autocorrelation  function 
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Then,, 


(B71a) 


(B71b). 


The  values  a ^  y  aa^  and  p  derived  from  Eq.  (B65)  ,  (B71a)  and 

(B71b)  can  be  used  in  the  joint  probability  statistical  analysis  described' 
previously. 
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APPENDIX  C 


NUMERICAL  RESULTS 


The  wing,  body  and  vertical  tail  loads  resulting  from  a  one  foot  per 
second  rms  gust,  A,  and  corresponding  number  of  crossings  per  second, 

N0>  are  tabulated'  oh  the  following  pages.  The  condition  numbers  cor¬ 
respond  with  those  shown  in  Tables  2  through  5  in  the  Analysis  Conditions 
section. 
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VERTICAL  ANALYSIS 
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